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Abstract 


The theory of Massera and Schaffer relating the existence 
of unique almost periodic solutions of an inhomogeneous linear 
equation to an exponential dichotomy for the homogeneous 
equation has been completely extended to discretizations by 
a strongly stable difference scheme. In addition it has been 
shown that the almost periodic sequence solution will converge 
to the differential equation solution at a rate O(k^) where p 
is the accuracy of the scheme, uniformly in t, if the coeffi- 
cients are sufficiently smooth. 

The preceding theory has also been applied to a class 
of exponentially stable partial differential equations to 
which one can apply the Hille-Yoshida Theorem. It is possible 
to prove the existence of unique almost periodic solutions of 
the inhomogeneous equation which can be approximated by 
almost periodic sequences which are the solutions to appropri- 
ate discretizations. Two methods of discretizations are 
discussed; the strongly stable scheme described above and 
the Lax-Wendroff scheme. 
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Introduction 


This work extends certain facets of the theory of 
Massera and Schaffer [1] (in the future referred to as M&S) 
relating properties of the solution of an inhomogeneous 
ordinary differential equation (ODE) to properties of the 
solutions of the homogeneous system. In Part I we extend 
this theory to difference approximations of the ODE. In 
Part II we use the ODE theory, exclusively in the stable 
case, together with the Hille-Yoshida Theorem to obtain 
results for partial differential equations. 

The study of admissibility theory and dichotomy theory 
for difference equations was first done by Coffman and 
Schaffer [2]. This work differs from their work in that we 
are concerned with the preservation of , admissibility and 
dichotomy properties by a difference equation which is used 
to approximate an ODE system which has certain of these 
properties. 

First let us review the theory of M&S. Consider the 
equations 

(0.1) y = A(t)y 

(0.2) y = A(t)y + f (t) 

Here the independent variable t ranges over the whole real 

line and for each t the vector y(t) lies in a Banach 

space E (which may be infinite dimensional) . A(t) , for 
each fixed t, belongs to the space of bounded operators on E 


- 1 - 



which we will denote by L(E,E). Let B and D be two 
Banach spaces of a function from the real line into E. 

We say the pair (B,D) is admissible for (0.2) if 

a) Vf ^ B 3a unique solution y (to (0.2)) e d, and 

(0.3) 

b) II yd D < Kl f I B . 

This definition is more restrictive than that of M&S. We are 
not going to deal with all the subtleties of their theory, 
but just with certain important parts. It is pointed out, 
however, that (0.3b) is actually superfluous (see M&S, 

Chapter 5) but we include it in the definition of admissi- 
bility for simplicity. 

Observe that (0.3b) states that if we write the y 
given in (0.3a) as 

(0-4) y = Ctf) 

then C, which is obviously a linear operator from B into D, 
is also bounded. We will be concerned almost exclusively 
with the case that D = B = A (E) where A is the space of 
almost periodic (AP) functions with range in E. The range 
in this notation will generally be omitted and we will simply 
write A, etc. , if no confusion might arise* At times we 

will also be concerned with the cases B = D = L or 

00 

6=0=0^ which we define as the space of bounded continu- 
ous functions with range in E. Observe that A and C m 
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are closed subspaces of and we have the inclusions 
(0.5) A C Cco C Lco 

A ( t) will always be an AP or operator function 
unless stated otherwise; 

Associated with the concept of admissibility we define 
the concept of an exponential dichotomy for (0.1). Specifi- 
cally (0.1) has an exponential dichotomy if there exist 
projections P 1# P 2 = I - such that if Y{t) is the 

fundamental operator solution to (0.1) (Y(0) = I, see M&S, 

Chapter 3) , the following estimates hold for some a > 0. 


a) 

1! Y(t)P 1 Y 

1 (s)» 

< Ke' 0(t ' s> , 

t < S 

b) 

II Y(t)P 2 Y~ 

1 (s) II 

<Ke- a(5 - t > , 

S > t 


Note that this is not the general definition of an 
exponential dichotomy given in M&S (Chapter 8) but it is 
equivalent to their definition when E is finite dimensional. 
In the case that E is infinite dimensional we will in 
general only be concerned with the stable case; that is, 

P 1 = 1 ' P 2 = 

Observe that P^ is merely the projection on E + , the 
space of initial data of solutions to (0.1) which are bounded 
for t € [0,<»). To see this note that if y Q were in the 
range of P 2 and if the corresponding solution y(s) were 
bounded for s >_ 0 then, by setting t = 0 in (0.6b), we 
obtain 
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(0.7) 


y 0 = Y _1 (s)y(s) = P 2 Y -1 (s)y(s) 

and Yq = 0 follows from letting s Similarly P 2 is 

just the projection on E_ , the space of initial data which 
are bounded for t e (-®,0]. An exponential dichotomy merely 
states that E = E + © E_ and that the solutions, in addition 
to being bounded, decay exponentially. 

Now it is a fundamental result of M&S that if A(t) e A 
then (A, A) is admissible for (0.2) iff (0.1) has an expon- 
ential dichotomy. The preceding statement is also true if 
A is replaced everywhere by L^. We prove here the easy part 
of the result; namely if we have an exponential dichotomy, 
then we have admissibility. If f e define 


00 


(0.8) 

y ( t) = c ( f ) = f G(t,s) 

J 

f(s) ds 


where 

— co 

G(t , s) = YttjPjY -1 

(s) , 

t > s 

(0.9) 

= -Y(t)P 2 Y -1 

(s) , 

s > t 


Observe that the estimate 

(0.10) HG(t,s)H <Ke“ a l t " s ' 

holds by the definition of an exponential dichotomy. 

Using (0.10) we see that the integral in (0.8) exists, 
that y satisfies (0.2), and that y is the unique L ffl 
solution to (0.2) and in fact 
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( 0 . 11 ) 


|( y |] <2K Bf || 

J 00 00 

So that C, defined in (0.8), is a bounded operator from 
(E) L^E). To see that C: A -*• A we merely let T 

be a common e-almost period for f and A. Then w(t) = y(t+T) 

- y (t) is the unique solution to (0.2) with 

inhomogeneous term 

(0.12) [A(t+T)-A(t) ]y(t+T) + f(t+T) - f(t) = 0(e) 

We then use (0.11) to obtain 
(0.13) i (y(t+T) - y(t))l„ = 0(e) . 

The proof of the converse is more difficult and can be 
found in M&S, Chapter 10. A simplified proof valid only in 
the finite dimensional case is given in the Appendix. 

In Section 1 we will define almost periodic sequences, 
which will be the type of solution we will be searching for. 
Section 2 is the most important of this work. Here, after 
discussing the properties of the strongly stable difference 
schemes we will be using, we will introduce a transformation 
which will separate out the roots of the scheme inside the 
unit circle and permit us to work with the Euler 1-step scheme. 
This technique was originally developed by Engquist [3] 
although the author was not aware of his work when the 
formulation given in Section 2 was developed. Engquist' s 
results will be discussed more thoroughly at the end of 
Section 3. 
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In Section 3 we settle the question of admissibility 
when the homogeneous system is exponentially stable. Using 
an inequality which is an exact discrete analogue to the 
Gronwall inequality, we will show that the homogeneous 
difference equation is also exponentially stable. From 
there it will be a simple matter to obtain, for any suffi- 
ciently small time step k, the existence of an almost 
periodic sequence as a solution to the inhomogeneous 
difference equation. Furthermore the sequences converge 
uniformly to the unique AP solution to the ODE with a 
uniform error 0(k p ) where p is the order of accuracy of 
the scheme. We will also show that this solution can in 
fact be calculated, i.e. it is stable under roundoff 
errors and errors in initial data. 

In Sections 4 and 5 we deal with the case that the 
homogeneous system has a general exponential dichotomy. 

We will show that the corresponding inhomogeneous difference 
equation also has an exponential dichotomy. This is only 
of theoretical interest as the solution will no longer be 
stable under roundoff errors or errors in initial data. 

In Section 6 we will deal with some miscellaneous topics, 
especially the convergence of the mean value of the AP 
sequence to the mean value of the AP solution, and also the 
weakly nonlinear case. 

In Part II we extend this theory to a simple class of 
partial differential equations which can be written as an 
evolution equation 
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(0.14) y = (B-5 ]y + f 

where B is an unbounded operator which satisfies the 
conditions of the Hille-Yoshida Theorem and 6 (t) is an 
AP function such that the homogeneous system is exponentially 
stable. We can obtain a unique AP solution to (0.14) by 
using the formula which would be valid if B were bounded 
and then showing that under mild restrictions on f the 
resultant function does in fact satisfy (0.14). 

In Section 8 we introduce a family of bounded operators 
which are spatial discretizations to B. We construct 
functions y h which are the unique AP solutions to 

(0.15) y h = [B h -6]y h + f • 

We will give conditions to insure that 

(o.i6) B <y-y h ) 11 m - °(h^) , 

where j is the order of the approximation of B h to B. 

In Section 9 we apply the theory of Part I to (0.15) to 
obtain an AP sequence y n;h which approximates y h - 
We have however the unfortunate restriction 

(0.17) JT5 = 0(1 > 

h 

where U B. II = 0{l/h m ). A more favorable result 
h 

(0.18) — 0(1) 

h 

is obtained in Section 10 when using the Lax-Wendroff scheme. 
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In Section 11 we will consider the extension of these 
results in the case that the operator B is perturbed by 
some bounded AP perturbation D(t). 

Finally in the Appendix we will give a proof of the 
basic ODE theorem that admissibility is equivalent to the 
existence of an exponential dichotomy. The proof is valid 
only in the finite dimensional case but is simpler than the 
proof given in M&S and is also simpler than a finite 
dimensional proof to be found in Coppel [4] . 
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1 . Almost Periodic Sequences 


Our first task is to introduce the discrete analogue 
of an AP function, Following Ccrduneanu [5] we define an 
almost periodic (AP) sequence a^ with range in E as 
follows : 

a n is AP iff given e > 0 there exists a length L(e) 
(a positive integer) such that in any sequence of L 
consecutive integers there exists an N such that 


( 1 . 1 ) 


fl a 


n+N 


- a II 

n 00 


< e 


where the sup in (1.1)' is taken over n. As shown in 

Corduneanu (page 45) this is equivalent to normality i.e. 

given any sequence of integers the sequence 

b . = a will have a uniformly convergent subsequence . 

n;i n+N^ J 

Although the proof given in Corduneanu is stated for scalar 
valued sequences this proof is obviously valid if the range 
is any Banach space E. 

If we define the space l£(E) as the Banach space of 
bounded sequences with range in E then the AP sequences A n (E) 
form a closed subspace of L^(E). As usual the argument E 
will be omitted when no confusion can arise. 

We point out that normality can be used, exactly as in 
the continuous case, to show that for any finite set of AP 

1 r * « , 

sequences a n '** ,,a n ^ange in possibly different spaces 

E and for any e > 0, there is always a length L(e) 
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such that in any interval of length L we can find a common 
e-almost period. 

Finally we note that if f(t) is an AP function then 
the sequence f = f(nk) is an AP sequence for any real k. 
The converse is also true as shown by Corduneanu (page 47) 
but we shall not use that. 
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2., Properties of the Difference Scheme 


We consider linear 5,-step, strongly stable schemes 
described as follows. 

(2,1) j l Q “j y n+j = k 3 j y n+j " k j| 0 B j [A n+j y n+j +f n+j } ' 


Here we assume that we are discretizing the ODE 

(2.2) y = A (t ) y + f(t) 

and y = y(nk), A = A(nk), f = f(nk) where k is the 
n n n 

time step. 

Associated with (2.1) we have the polynomials 

St, i . ' 

(2.3) p(x) = . I a.x- 1 , a (x) - l g.x 3 . 

j=0 j=0 

It is well known (see Dahlquist [6] or Heinrici [73) that 
consistency implies that x = 1 is a simple root of p(x) = 0 
and that 

(2.4) p' (1) = o(l) = 1 

where we have normalized the coefficients so that the common 
value in (2.4) is one. 

Strong stability of the scheme is achieved by restricting 
the size of the other 5,-1 roots of p(x) = 0. Specifically, 
if we number these roots x y , u = 1 then we require 
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that there be a positive number 0 < 1, such that 

(2.5) - 1 , | x u | < 9 < 1 , u* 2, 

Finally, we make two further assumptions. 

(2.6) distinct , u = , 

(2.7) / 0 , u — 2 . 

(2.7) is necessary because to get an AP solution, the 
difference equation should be solved backwards and forwards. 
This assumption can be removed in the important case that 
the homogeneous ODE system is exponentially stable. Condition 
(2.6) can be removed in all cases and is included here only 
to simplify the following proofs. The removal of these 
conditions will be discussed in Section 6. 

Now in working with a multistep scheme, the standard 
procedure is to convert it into a one-step scheme. To do 
this we define the space E^ = Exex...xe. (We will usually 
write vectors in E^ in column vector form. ) We give E^ the 
norm inherited from this definition, namely if w € and 



then 

(2.9) 8 wii ^ = max By^ilg . 

i 

Here we have explicitly indicated the E norm in (2.9). 
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We now consider the discretization of the linear 


inhomogeneous system (2.2). If = 


^ n +&- 1 

i ) 


then we ge t 


( 2 . 10 ) 


w . = C w + kf . 
n+l n n n 


Here if we define z = k A , 

n n 


(2.1!) C n - C(z n# ... ,z n+a ) - 


* a * I-e jl z n+A* (6 S,-l z n+J>-l~ a £-l I) 


•• (a l I ” B J 2 n+it ) ' 1(e 0 2 n" a 0 I) 


I 0 ... 0 

0 1 . : 

* 

♦ 

* 

: i 

o 


0 ... 0 


0 10 



f 


n 


^l 1 ~ e £ z n+A ) 6 j f n+j 

0 

* 

* < 

m 1 

o 


Observe that C e L(E„ / E /> ) and will be just an mixmJl 
n x. x 

matrix in the case that E is an m-dimensional space (in which 
case it is called the companion matrix). Note also that since 
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A is uniformly bounded in t, the implicit term (ot^l-B^z^^) 
can be inverted, for k sufficiently small, uniformly in n. 
We finally point out that if A and f are AP (L^) then 
and are AP (L^) and also that C given in (2.11) is a 
smooth function (in the Frechet sense) of its &+1 arguments. 

Now with A a constant and z = kA € L(E,E) we consider 
for small B zB the homogeneous difference scheme 

(2.12) w n+1 = U(z)w n 

U(z) = (ct £ i-3 A z)‘ 1 (B Jl _ 1 z-a £ _ 1 I) ... (a £ I-3 Jl z)“ 1 (B Q 2-a 0 I) 

10 .. . 0 

0 * . * . 

1 0 , 

« » • 

• • 

0 

0 01 0 


U(z) is a mapping, defined for small II zll , from L (E,E) into 
L(E^,E^). Observe the following properties of U(z), 




Now (a) and (b) should be clear. To see (c) , note that the 


spaces E are closed and also that any two have only zero 

?1 

in common. Let w = I and suppose w has an expansion 


(2.14) 


x £ 

l 

w = l 

U— 1 


fx^y 

u . u 


x y 
ir u 


l 


u 


Then if we define the vector w to be 

4 

(2.15) w = V w 


t y l 


we can write 


V 


- x. 


£-1 


x 


l-l 

i 


X 1 ’ x £ 

1 ... 1 

i.e. V is just the Vandermonde matrix associated with the 

distinct numbers (2.13c) now follows directly 

from the invertibility of V. 

The representations (2.14) and (2.15) enable us to 

define an equivalent norm on the space E^. Specifically if 
w is expressed as in (2.14) we define 


(2.16) il w tl u = max 1 y II E . 

u 

u o',.. 

That K II is equivalent to II II follows immediately from 

the representation (2.15) and the invertibility of V. Note 

that the equivalence of these norms implies the equivalence 

of the operator norms they induce on L(E^,E^). Of course 
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this paragraph is superfluous in the case that E is finite 
dimensional. 

For future use we point out that (2.13b,c) imply 

that u" 1 (z) exists for small II zll and is smooth in z. 

Now define the space E^ = © £ E u « We have E^= E.©E^ 

u> 1 

and we note that is canonically isomorphic to E. At times 
we will identify E ^ with E but this should not cause any 
confusion. With respect to the decomposition E = E^ E 1 we 
see that U(0) is in block diagonal form; symbolically, 

f D{0) 0 'i 

(2.17) U(0) = 

^ 0 B (0 ) 

here D(0) : E^ + E^ and is the identity, while B(0); 
and H B ( 0 ) 11 U < 0. This is the operator norm induced on 
L(E^,E^) by the II II U norm on E^ and follows from the fact 
that B(0) is just multiplication by on the space E^. 

We can now state the fundamental theorem of this section. 

Theorem 1 . For small II zll there exists an operator 
T (z) : L(E,E) -> LfE^E^) such that 

(a) T(0) = I (Identity on E^) 

(b) T(z), T 1 (z) are smooth in z 

_ i 

(c) L ( z) = T (z) U(z) T (z) is in block diagonal form 

^ 1 

with respect to the decomposition © E . 

(d) Writing L(z) symbolically as 

fD(z) 0 'i 

L(z) = 

1 0 B (z) ; 
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then B ( z) ] = B ( 0) (from (2.17)) and llB(z)ll u < 6 while 

Z—“ u 

D (z) has an expansion 

2 

(2.18) D(z) =1+ z + 0 (I zll ) 

where is identified with E. 


Theorem 1 is basically trivial and the proof involves 
familiar arguments. First we note that U(z) is smoothly 
invertible for small I zll . (The restriction " for small H zll " 


will not be stated explicitly in the future.) 

Next we observe that (wl - U(z)) exists for w in a 
small annulus around the circle |w | =0 and this holds 
uniformly in z. This follows from (2.13b) and the geometric 
series . 

Now define the projections 


(2.19) P(z) = j 


(wl-u(z)) 1 ,dw , Q(z) = I - P (z) . 


, W 


| = 6 


That P and Q are projections is a familiar result which 


follows from the resolvent identity, 


(w x I-U(z)) 1 - (w 2 I-U(z)) 3 


( 2 . 20 ) 


= (w 2 -w x ) (w 1 I-U(z)) 1 (w 2 I- U (z) ) 1 , 


2 . . 
calculating P by integrating around two slightly different 

circles, and interchanging the order of integration. It 

follows from the construction that P(z) and Q(z) commute 

with U(z) and that they are smooth functions of z. 
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It is clear by applying P{0) to an arbitrary vector 
w £ E £ and using the expansion given in {2.15) that P(0) 
is exactly the projection onto E^ along and Q(0) is the 
projection onto E^ along e\ 

Now define 

(2.21) T ( z) = P (z) P (0) + Q ( z ) Q ( 0) . 

Clearly T(0) = I and T(z) is smooth whence T ^(z) exists 
and is also smooth. Note that this holds in either of the 
norms U II U or II II £ on L(E^,E^). Let 

(2.22) L ( z) = T~ 1 (z)U(z)T(z) . 

We claim L is in block diagonal form with respect to 
the decomposition E^ - E^ © E^ . This follows directly 
from the fact that P(z) and Q(z) commute with U(z). In 
fact let x s E^ ; then 

(2.23) U(z)T(z)x = U(z)Q(z)x = Q(z)U(z)x = w (say). 

Now if y = T ^(z)w and y = y^ + y^ with y^ e 
1 1 

and y e E we have 

(2.24) w = Q(z)U(z)x » T(z)y = Q(z)y^ + P(z)y^ 
whence 

(2.25) P(z)y 1 = T(z)y^ = 0 

and so y"*" = 0 by the invertibility of T(z). Thus L(z) 
maps -> and a similar argument shows L(z) : E^. 
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If we write {using the notation of (2.17) ) L(z) as 

D ( z) 0 D ( z) : ® i 

(2.26) L ( z) — 1 i 

0 B ( z) B (z) : E E 

then we see that D(z) and B(z) are smooth in z; D(0) = I 
and Ab( 0)!I U < 0. It follows immediately by continuity 
that 

(2.27) Ob(z)!I u < 9 


for small II z II . Now (2.27) will imply that given z^#...#z^ 
With II z J small then 

(2.28) » 1 T B(z )ll £ < K e N 

i=l 1 

for a constant K independent of the z^ and N. 

It only, remains to study the first order structure of 
D ( z) in order to obtain (2.18). Expanding D(z) about z = 0 
we can write 

(2.29) D ( z) = I + C { z) + 0 { II z II 2 ) 


where C is the Frechet derivative of D(z) at z = 0 and is 
a bounded linear map from L(E,E) into L(E-^ / E 1 ) which we 
identify with L(E,E). It is necessary to show that C is 
the identity, and as one might expect this will follow very 
easily from consistency. 

Let A be an arbitrary element in L (E , E ) and let z = kA, 
A* = C (A) . For y € e and identifying E-^ with E, we 
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calculate the following limit 

(2.30) lim D(kA) n y = (I + kA ’ +0(k 2 )) n y 
k ->-0 

nk= t 

for any fixed t > 0 . 

tA ' 

As one would expect this limit is e y because 

the difference scheme expressed in (2.30) is consistent 
* 

with the ODE y = A'y. 

More precisely since 

(2.31) I D (kA) II < 1 + kR 

for small k f where R is some fixed constant, it follows 
that 

(2.32) JJ D (kA) N !f < e RT , Nk < T, 0 < t < T . 


Let x_ 
n 

(2.33) 



i.e. x = D(kA) n y. Now if x = e A ^ n y then satisfies 

n n n 

the perturbed difference equation 


(2.34) x , = [D (kA) + 0(k 2 )]x„ 

n + 1 n 

2 

where the term 0(k ) is uniform for nk 5 T. Equation (2.34) 
follows directly from the fact that x(t) satisfies the ODE 


(2.35) 


x = A* x . 
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Letting w n = x r - x n we can sh 9 w eas i!y enough that 

(2.36) II w II = 0 (k) , nk < T . 

n — 

In fact from (2.33) and (2.34) we obtain 

7 n i 

(2.37) w n = k l D(kA) 3 0(l)w .. 

n j-1 J 

If we use (2.32) to bound the powers of D(kA) and thus H x n ll , 
we then obtain (2.36) and letting k -► 0 we obtain 

n tA ' 

(2.38) lim D(kA) y = e y 

n-*-°° 

nk=t 

Recalling the identification of with E we define the 
Y 

vector y = ! whence y e E-.. 

y 

Thus (2.38) can be expressed 

(2.39) lim L(kA) n y = 

k-»-0 

nk= t 

Now we are going to show that A' = A by calculating this 
limit directly from the definition of L(kA). In fact since 
U U ( 0 ) II U = 1 it follows that »U(kA) n H u is uniformly bounded 
for nk _< T (and of course this is also true if we replace 
II 8 u by II I! 5 '). Then writing 

(2.40) L(kA) n y = T _1 {kA) U(kA) n T (kA) y , 

and letting k •+ 0 such that nk ® t, using the continuity 
of T ( z) and the fact that T(0) = I together with the 
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uniform boundedness of the powers U(kA) n , 

(2.41) lim L(kA) n y = lim U(kA) n y 

k+0 k-*0 

nk =t nk=t 

the last equality follows from the fact that a consistent 
and stable difference scheme is convergent. 

Since A was arbitrary, (2.41) together with (2.39) shows 
that A = a', which proves that the map C is the identity 
and this completes the proof of Theorem 1. 
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3. The Exponentially Stable Case 


We now assume that the homogeneous system (0.1) is 
exponentially stable. Thus, if Y(t) is the f undamental 
operator solution to (0.1) (Y(Q) = I) , then the estimate 

(3.1) llY (t)Y _1 (s) II < Ke 0<(t " s) , t >_ s , 

will hold with positive constants K and a. In the sequel 
we will use K as a generic positive constant so it will 
appear in contexts other than (3.1) . 

We wish to prove a similar estimate for the homogeneous 
difference equation (see (2.10), (2.11)) 


(3.2) 


w i = C w . 
n+1 n n 


Wow C is invertible for all n by (2.7) . Then (3.2) 
n 

has a fundamental solution V? n (Wg = I) which is simply 


■ . n-1 

W n - FT C i ' 


(3.3) 


i=0 

-n 


W 


-n 


i=-l 


C r _1 ' 


n > 0 , 


n > 0 . 


It can easily be verified that the unique solution 
to (3.2), w n , given initial data w^ is simply 


(3.4) 


w = W W.^-w. 
n n j j 


Note that for n > j, w n w ^ 


-1 


is simply 


-2 3- 



(3.5) w w " 1 - T~r C. , n > j . 

J i=j 

In the stable case we will only use W n W_.^ f° r n il j 
and we can see that the invertibility of is not required 

in this case and assumption (2.7) is therefore not required 
here . 

We intend to show that for small k (this will not be 
stated explicitly in the future) (3.1) implies 

, -a.k(n-j) 

(3.6) lw n w j * £ K e , n > j , 

where K is used as a generic constant and can be taken 
independently of k, while == a + 0(k). Thus since a > 0, 

can also be taken independent of k for small k. 

To attain (3.6) we suppose A(t) is C^. All derivatives 
that are assumed will always be required to be AP or L^. This 
will not be stated explicitly in the future. The smoothness 
condition on A will be removed in Section 6. 

Expanding A n +j about A^ we can write (3.2) as 

(3.7) w i = U (kA )w + k 2 0 (1) w . 

n+i n n n 

Here the 0(1) term is uniform in n and the operator U is 

given by (2.12). We are going to show that perturbations of 
2 

the type k 0(1) preserve exponential stability. We thus 
consider the unperturbed system 

(3.8) 5 n+1 = U(KA n >W n 
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If we let v n = T(kA n )w n , where the operator T is 

introduced from Theorem 1, and note that by the smoothness 

2 

of T and A we have T(kA„ . , ) = T (kA ) + k 0(1), we see that 

n+ i n 1 

we should consider the difference equation 


(3.9) 


v 


n+1 


L (kA ) v 
n n 


where the operator L has been studied in Theorem X. 

If V denotes the fundamental solution to (3.9) then 
n 

since L is block diagonal we see that V n itself is also 
block diagonal. If we write V n in the form 


(3.10) 


V 


n 



n 


l 0 v " J 


we can see that the second part, the contribution from the 
roots inside the unit circle, causes no difficulty. In fact 
for n > j we have 


o --1 n-1 

vV - TT B OcA. ) 
n 3 i=j 1 

(3.11) 

ii v 2 v2~ 1 ii < k e n-j , 
n 3 


where we have used (2.28) and K as a generic constant. 
Since 0 is a fixed number less than 1 we see immediately 
that (3.6) holds for where can in fact be taken 

equal to a for small k. 

We must now deal with V^. Neglecting perturbations 
of k 2 0(l) and identifying the space E 1 with E we 
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consider the difference equation 


(3.12) x n+1 = (I + kAn)x n 

i.e. we have reduced the problem to the study of the Euler 
1-step scheme. 

We will compare X n , the fundamental solution to (3.12), 

with Y (= Y(nk)) where Y is the fundamental solution to 
n 

the homogeneous equation (0.1). 

Let H(t,s) = Y (t) Y ^(s). As a function of t, 

H satisfies the equations 

H = A ( t) H , 

(3.13) H = A 2 (t)H + A ( t) H , 

H ( t, t) = I 

Integrating (3.13) we obtain 

1 

H((n+l)k,nk) = I + k | H(nk+0k,nk) d0 

(3.14) 1 1 0 

= I + kA(nk) + k 2 j j 0 d0 H(nk + <t>6k, nk) . 

0 0 

Now since A is bounded it is a standard result (M&S, 

Theorem 3.1 C) , that 

(3.15) «H(t,s)ll < K X (R) for |t-s| < R , 

where R is any positive number. This does not require 
exponential stability and in fact is an immediate consequence 
of Gronwall's inequality. It now follows from the second 
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2 

equation in (3.13) that the coefficient of k in (3.14) is 
bounded uniformly in n and k. This implies that Y fi is 
the fundamental solution operator to an equation 

(3.16) Y n+1 = [I+kD n (k)]Y n = [I+kA n + k 2 0(l)]Y n 

where the 0(1) term is uniform in n. This means that we 
can regard the x^ equation (3.12) as a perturbation of the 
Y^ equation. We point out that the reduction we have obtained 
does not use exponential stability and so is valid in the case 
that (0.1) has a general exponential dichotomy. If we use 
the stability, however, we see that Y n satisfies the 
estimate 

(3.17) II Y. yT 1 H < k e -ak(n "^ , n > j . 

n 2 — 

We must now prove the proposition that exponential stability 

2 

is preserved under k perturbations. 

Theorem 2 . Consider two difference equations 

(3.18) y n+1 = R(n,k)y n 

(3.19) x n+1 = R(n,k)x n + kS(n,k)x n 

defined for small k. Suppose that (3.18) is exponentially 
stable, that is there exists constants K and a independent 
of k such that 

(3.20) 8 Y^yT 1 ' 1 < K e ~ ak (n- j ) ^ n > j , 
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where Y is the fundamental solution to (3.18). We point 
n 

out that we do not require R(n,k) to be invertible since 

we only study the solution Y^Y^ for n j. In that case 

Y is not invertible and we must use (3.5) in place of Y Y ^ , 
n n j 

but we will retain this notation for simplicity. 

Under this hypothesis there exists an such that if 

for small k. 


(3.21) 


H S (n ,k) II £ e ± £q 


where the sup in (3.21) is over n, then (3.19) is exponentially 

stable and if X is the fundamental solution to (3.19) (the 

n 

dependence on k has been suppressed) the estimate 

, -a.k(n-j) 

(3.22) II X n Xj II £ f^e 


where can be taken independent of k and = a + 0(e). 

Before proving this we note that the case where S(n,k) 

= 0(k) is automatically covered. The more general formula- 
tion will be used in Section 6 and more importantly in Part II. 
Note also that it is a discrete analogue to the Gronwall 
inequality and the proof is in fact immediately suggested by 
the proof of the Gronwall inequality. This theorem is equivalent 

to a lemma of Engquist itself based on a theorem of Strang [8], 
but the proof given here is simpler than Strang's proof and 
much more suggestive of the Gronwall inequality. 

To prove Theorem 2 we first note that for n > j we have 


(3.23) X n xT 1 = 


Y yT 1 + k 
n 3 


l 


£=j + l 


~ — 1 
Y Y ■ 
n 3 


S 


£-1 
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This analogue of the variation of constants formula can be 
verified immediately. Using (3.20) we get 


(3.24) II X X . ~ 1 II 
n 3 


n -1 


< K e 


-cck(n-j) + kK ak r & -ak(n-L) x' 1 ! . 

lij * ’ 


Here K is the constant of (3.20). If we redefine K as 

K e 01 ^ (for small k) and define v . = e a ^ n lx X.^l we 

n, 3 n j 


can write 


n -1 

(3.25) v nij <K + k Ke l v tf) = < n>j 

X,— _| 


From (3.25) we obtain 


(3.26) 

whence 


v 


n,3 


s .-s . . 

= ,D /.3 < s 

kKe - n-1, j 


(3.27) . n)j ' < (1+eKk)" - -^ s j (j 5 *i S EKk(n ' j) 


where K. is defined so that s. . <_ K, independent of k 
-L 3 f J -L 

and j (from (3.25) specialized to n = j+1 and the fact that 

v . . = 1) . 

JfJ 

Inequality (3.27) together with (3.25) and the defini- 
tion of v . yields (3.22) immediately and thus completes 
n I 3 

the proof of Theorem 2. 

Returning to our specific case, it should be clear that 
after several applications of Theorem 2 we can prove 
exponential stability for v R defined as T(kAn)w n (see (3.8) 
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and the following paragraph) . Since in terms of fundamental 
solution operators we have 

_ _-l -1 

(3.2 8 ) w W. = T (kAn)V V. T(kA.) 

n 3 n 3 3 

we obtain exponential stability for (3.8) and another applica 
tion of Theorem 2 yields exponential stability for the full 
homogeneous system (3.2). It is a simple matter to go from 
exponential stability to admissibility for the inhomogeneous 
difference equation 


(3.29) 


w , . = C w + kg 
n +1 n n n 


We first define admissibility in the obvious way. 

The pair (L^,L^) will be admissible for (3.29), for small 
k, iff for any sequence g^ in i/J (3.29) has a unique solu- 
tion w n in and this assignment is a bounded mapping, 

i.e. there exists a K (independent of k) such that 

(3.30) llwll < Kl! g I . 

n °° — n 03 


An entirely analogous definition holds for the admissibility 
of the pair (A n ,A n ) . As in the ODE case we will show (L^rL^) 
is admissible and then show that the solution is AP if the 
coefficients are. 

To show (L n ,L n ) admissibility we merely write the solution 


(3. 31) 


ft T_ 

w = k T W W . g . . . 
n > „ n 3 3-1 


l=- c 


Exponential stability implies the convergence of this series. 
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which can be immediately verified to be a solution of (3.29). 
To show w n £ we estimate, using (3.6) 


(3.32) Itw I < II g II Kk V e“ ajk 
n 00 — “> 


g II 
y n 00 


K 


1- e 


■ak 


< II g « K 
— 3 n “ 


(using x/(l-e x ) ■+ 1 as x -*■ 0) and this yields (3.30) after 
a redefinition of K. 

To show uniqueness is trivial, if were a bounded 

solution to the homogeneous equation (3.2) we would have 
for n >_ j , 

(3. 33) w = W wT 1 w . 

n n j D 

and w n = 0 follows immediately upon letting j -*• 

It remains to show that (A n ,A n ) is admissible. In 

fact if N is a common e-almost period for and g n , 

we would have w , - w to be the unique L solution 

n+N n 00 

to an inhomogeneous equation with inhomogeneous term 

(3 ' 34) k fC n+N~ C n^ w n+N + g n+N ” 9 n 

and for fixed k > 0 the almost periodicity of w^ is a 
consequence of (3.30). 

We can remove the factor 1/k in the denominator by using 
the expression (2.11) for C > 

(3.35) C n = C(kA n ,. . . ,kA n+k ) , 

where C is a smooth function of its arguments. Expanding 
C n+N _C n us: *- n 9 (3.35) we see that if N is a common e-almost 
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period for A and f , then N will be an 0(e) -almost period 
n n 

of w , with the bound in the 0(e) term independent of k. 
n 

Note that these two paragraphs follow only from (l£,L^) 
admissibility of (3.29) and do not depend on the stability 
of (0.1). 

The convergence of w n to the solution y n also follows 
very easily. If y is the unique AP solution to (0.2) we 
have 

<3. 36) Z a jyn+ .-k [ 6 j y n+j 

1 

= k[ £ a . j j (y (nk+6k j) - y(nk)) d6 

3 0 

- I 6 . [y (nk+kj) - y(nk)]] = k g(k,nk) . 

j J 

(The summation is from 0 to & in the above.) Here g(k,t) is 
AP in t for each k and 


(3. 37) 


g(k,t)» = o(l) , 


k -► 0 ; 


(3.37) follows from the uniform continuity of AP functions on 


the whole real axis. 


n+£-l 


then (3.36) implies that w n is the 


Now if w =i . 

n ' y 

unique AP solution n to the same equation as w n except for 


an error. 


(3.38) 


f g(k,nk) 
0 
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and (3.37) together with (3.30) yield 


(3.39) II w -wl = o(l) , k -► 0 . 

n n °° 

We can obtain finer convergence results by imposing some 
smoothness conditions on A(t) and f(t). In fact if the scheme 
has order of accuracy p, and A and f have p+1 derivatives, 
then y also has p+1 derivatives. It is shown in Henrici 
(p. 247) that 

(3.40) II l ajy n+j -k l 6jy n+j * < k p+1 Glly P+1 ll co , 

where the constant G depends only on the scheme. (3.40) 
together with (3.30) yield 

(3.41) lw n -w n y OT = 0(k p ) . 

We point out that these convergence arguments are valid when- 
ever we have admissibility and do not require stability. 

Before leaving the stable case we would like to discuss 
problems relating to the computability of the solution w n - 
Consider first errors in initial data. Suppose we 
solve the exact difference equation (3.29) but use as initial 
data w*(k) = w Q (k) + e(k) where w Q (k) is the exact initial 
data for the AP sequence solution w n and e (k) is bounded for 
small k, i.e. le(k)l < e Q . Now the solution w* , which we 
solve for, will be 

(3.42) w* = w n + W n e(k) , 
and we have 
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ft — rv n 

(3.43) * w n~ w n* " °( e °) t n 0 . 

Now in general we cannot expect e(k) + 0 as k 0, 

This would certainly occur if we could obtain initial data 
consistent with the initial data of the solution y to the ODE. 
However 'this solution is known only by an improper integral 
involving not only the given forcing function f(s) but also 
the fundamental solution Y(t) to the homogeneous system 
which we could not expect to know explicitly unless A were 
a constant. Thus the most we can assert is that, using fixed 

ft 

initial data (say w^(k) =0), if we integrate over a sufficiently 
long interval we will get the solution w n up to some exponen- 
tially decaying error. 

Now let us consider round-off error. If we could 
postulate that we solved a perturbed equation 

* ft ft ft ~ r 4- 1 

(3.44) w n+l = C n w n + kf n + J • n t. 0, r > 0 , 

where the error is uniform in n then we could relate this 

ft 

solution to w by 
n 1 

ft ft ft ^ __ I 

(3.45) w = w + k l W w7 0(k r ) , n > 0 , 

n n , L . n j 

3=1 

and we can estimate the error simply, 

n 00 

(3.46) Ilk l W W7 1 0(k r )l _< k r Kk £ e“ akr =0(k r ) . 

j=l n 3 r=0 

Let us further examine the validity of (3.44). If we 
include rounding errors in forming the companion matrix (or 
operator) , together with the error of taking only a finite 
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expansion of the term (a^I ~ 6^kA n+ ^) ^ in the implicit 
case (see (2.11)), then we could assume that the linear 
part of (3.44) would be replaced by 

* * r+ 1 

(3.47) C = C + 0(k tTi ) . 

n n 

If we also assume an error 0 (k llw^ll) in forming the 
product C n w n we see we should postulate a system 

(3.48) w** = C w** + kg + 0 (k r+1 ) + 0 (k r+1 U w **ll ) 

n+1 n n n 11 


and in order to justify (3-44) we must show that the solution 
to (3-48) will be bounded for n 0 . 

We can write this solution as 

n 


** ** ri -1 

w n = W n w 0 + k .| x "n W j g j-l 


13.49) 


n 


„ _ i ~r -r * * 

+ k[ l W n Wj(0(k ) + 0(k •w ; j_ 1 l)] 


Now if we use the fact that 
(3.50) 


k l e akr = 0(1) 


r=0 


and define h = max 


** 


II w. 0 we can estimate h by 
n j=0 , . . . ,n 3 n 


(3.51) h n = 0(1) + k r h n 0(1) 

where the 0(1) terms are independent of n. This implies, 

, ^ ^ 

for small k, the boundedness of lw II and so justifies (3.44) 
The techniques used in Section 2 of block diagonalizating 
the companion matrix were first developed by Engquist 
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although I was not aware of his work when my formulation was 
developed. 

Engquist’s procedure may appear to be more complicated 

than the procedure presented here, but this is because he 

proves a more general result; namely that the companion 

matrix can be block diagonalized to within an error of 

p l 

arbitrary order (say 0(k )). (He also does not restrict 

himself to strongly stable schemes, but this imposes stability 

requirements on several additional homogeneous systems.) 

If one were to consider his procedure restricted to a block 

2 

diagonalization up to 0(k ) then, while the two formulations 
differ in the lines of approach, they are equally simple. 

Engquist also shows that the upper block will agree 
with the Taylor series expansion of (0.1), up to terms of 
order p, if p < is the accuracy of the scheme. The 

proof of the first order structure of the upper block (the 
term D(.z) in Section 2) is, I believe, somewhat simpler 
than his proof restricted to the first order term. 

I would like to point out that in the case we are 
considering (an exponentially stable ODE and a strongly stable 
scheme) the extra fineness of the block diagonalization being 
carried to order p+1 will only give the advantage of the 
homogeneous system having a stability exponent -a+0(k p ) 
instead of -a+0(k). In particular in Engquist's study of 
the uniform convergence (for t >_ 0) of the solutions of the 
homogeneous equation to the solution to the ODE (see his 
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Theorem 3, p. 24-27 and note that what he uses as a is 
what we call -a) one would still obtain uniform convergence 
of 0(k p ) (neglecting round-off and errors in initial data), 
however the decay exponent would be -a+0(k) instead of 
-a+0 (k* 3 ) . 
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4 . Admissibility in the Case of a 


General Exponential Dichotomy 


Here we are concerned with admissibility properties 
in the case that the homogeneous equation (0.1) has a 
general exponential dichotomy as described by (0.6) with 

0. The results will be valid in the infinite dimen- 
sional case, although in this case one must bear in mind 
that the situation described by equations (0.6) is not 
the most general form of an exponential dichotomy. 

First of all consider arbitrary homogeneous and 
inhomogeneous difference equations 


(4.1) 

(4.2) 


w n+l - C n <k,W n 


W„ , 1 = c (k)w + kg 

n+1 n n n 


Here C n is defined for k e (0,^] and is AP in n for 
fixed k, and is now assumed invertible so that a fundamental 
solution operator W n , as described by (3.3), exists and 
is invertible for all n. 

Associated with (4.2) we have variation of constants 
formulas in both the forward and backward directions. 


(4.5) 


(a) 


(b) 


n -1 

w ssWw^+k y W W. g. . 
n n 0 n j y 3 “l 


w n = Vo ' k jXiV'jS-l' 


n > 0 


n < 0 


as can be easily verified. 
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Now we will define an exponential dichotomy for (4.1) 
as follows. Suppose there exist projections P^(k), P 2 00 
- i - P^(k), such that the estimates 


(4.6) 


iw P.wT 1 !! 
n 1 3 


II W P„W . 1 ll 
n 2 j 


< K e-“ k(n -3> , 

< K e ‘“ k(j ‘ n) , 


n > j , 


j >. n • 


hold with K and a independent of k. The only difference 


between this and the ODE case is that the projections 


may depend on k. 


If we have an exponential dichotomy then (L^,!^) is 
admissible for (4.2). In fact the unique bounded solution 


w is 


(4.7) w (k) = k l <5g. , 

11 — — CO J J 


where 


G n = W P.W . 1 , n > j 

j n 1 3 — J 

(4.8) 

= -W P-W7 1 , j > n+1 . 

n 2 j J - 


Equation (4.7) is of course suggested by the ODE case and 
is easily verified using (4.6). The admissibility bound 
(3.30) can be derived exactly as in the stable case as can 
(A n ,A n ) admissibility. 

We next observe that (L^/L^) (or (A n ,A n )) admissibility 
is preserved under perturbations of the linear term. In fact 
if we had a system 
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(4.9) 


w n+l = C n (k)w n + kS n (k)w n + k9 n 
where, for small k, 

(4.10) l! S n (k)l oo < e 0 

and e Q is some number to be determined, then the unique l£ 
solution can be found by defining w n;0 E 0 and w n;i+i to be 
the unique solution to 


(4.11) w , . , = C w . . .. + kS w . + kg 
v'i.o.j.; n+1; l+l n n,i+l n n,x n 


(for simplicity the dependence of C n and on k has been 
suppressed). If £ () K 1 <1, where K 1 is the admissibility 
bound for (4.2), the contracting mapping principle establishes 
(L^,L^) admissibility with admissibility bound 


(4.12) 


K. 


K 2 (l-e Q K 1 ) 


Of course we will get AP solutions if D^+kS^ and g^ are AP . 

Now these two principles certainly settle the question 
of admissibility when the homogeneous system has a general 
exponential dichotomy. Referring to the reduction obtained 
in Section 3 we see that after a nonsingular change of 
dependent variable, = T(kA n )w n , v n satisfies the 

equation (see (3.8) f f. ) ; 


(4.13) v, 


n+1 


I+D (k) 0 


n 

0 B 


v + k 0 ( 1) v + kg 
n n 3 n 


n 


Here g n = T(kA n+1 )g n . B Y the boundedness of T and T 1 , 
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admissibility for (4.13) is equivalent to admissibility 
for (4.2). But (4.13) expresses the homogeneous term as 
a perturbation of the system 


( i+u u 
X . = [ n 

n+1 v o B 


with fundamental solution 




r Y 
f n 

0 1 




X = I 

n-1 

f 

n > 0 / 


n \ 

L 0 







i=l 



(4.15) 


y 

-n 

0 > 




X 


-n+1 n ! 

9 

3 

V 

o 


-n 

l 0 

FT b' 1 ! 

i=0 1 




and this has an exponential dichotomy with projections 
independent of k and in fact (after identifying with E) 


(4.16) 


P 1 = p x 8 p 




Here and P 2 are defined in (0.6) while P(0) is the 

projection onto E"^ along E, (see (2.19) ff) . 

We thus obtain admissibility in this case and we point 
out that the convergence proofs given in Section 3 did not 
require stability and remain valid here. 
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5. 


Existence of an Exponential Dichotomy for 


the Homogeneous Linear Equation 


We will prove here that admissibility for the inhomo- 
geneous difference equation 

f 5 - 1 * w n+l * C n"n + 

implies an exponential dichotomy for the homogeneous equa- 
tion 


(5.2) 


w n+ 1 C n W n 


The proof is very similar to a proof of the analogous property 
for the ODE case given in the Appendix. It is a strictly 
finite dimensional proof and it may then be simpler to think 
in terms of matrices rather than linear operators. We can 
then regard T(z) as a similarity transformation which follows 
an initial transformation putting U{0) into block diagonal 
form. The matrix L(z) is strictly in block diagonal form 
and it is no longer necessary for us to continually make the 
qualification "identifying with E." 

If we let v = T(kA )w , (5.1) is transformed into (see 

n n n 


(4.13)) 


(5.3) 


n+1 


I+kD 

n 

0 


0 

B 



+ k 2 0(l)v n + kg n . 


n 

(For simplicity we will write the inhomogeneous term in (5.13) 

as g rather than g . ) 
n 3 n 
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If we call the leading matrix in (5.3) R n (k), then the 
unperturbed linear system 


(5.4) 


x . . = R X 
n+1 n n 


has an exponential dichotomy with projections P-^ and P 2 / 
independent of k (see (4.15) and (4.16)) , and we must exhibit , 
an exponential dichotomy for the homogeneous version of (5. 3) . 


(5.5) 


n+1 


= R v 
n n 


+ k 0 (1) v 


n 


By admissibility (5.5) can have no solution. Let 

S + (k) and S_(k) be the subspaces of initial data which give 
rise to solutions bounded for positive and negative n respec- 
tively. We have S + ^ S_ = {0} (suppresssing the k dependence) 
If we let S be any complementary space we have 


(5.6) E a = S + © S_ © S 


with associated projections 

(5.7) 1 = Pi + P 2 + p 3 

Our first task is to show S = {0} i.e*. p 3 = °* 

Define L n to be the submanifold of of sequences 

with only finitely many nonzero components. Let g n e l^ 
and let v n be the corresponding unique L^ solution to 
(5.3). Then using (4.5a). we see that for large positive n 
we have 

(5.8) v n = V v 0 + k V j X g j-l ] 
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is bounded the vector in the brackets lies 


and since v, 


n 


in S + i.e. 


(a) P 2 V 0 = ' k .1 P 2 V J 1 9i-l 


(5.9) 


j=l 


(b) 


P 3 V 0 = 


“ ' k J, P 3 Vj 1 *j-l 
D=1 J 


Similarly looking at v r for n -*■ -« and using (4.5b) we 
obtain 


(5.10) 


(b) 


-1 


(a) P^o = k i J p ! V- g. x 


P v - k l P 3 v : 1 9j . 1 

j-: — co J J 


Now if we simply set 


g 0 - V 1 P 3 Z 


(5.11) 


g^ — 0 , i ¥■ 0 , 

where z is an arbitrary vector in we see immediately 
that P 3 = 0. We can also see that for any given g n e i£ 
the unique solution v n is given by 


(5.12) 


v = k l G n g . . 

n j^_<x> 3 


where Gj is given by (see (4.8)) 


(5.13) 


G 5 - V n P 1 V j X ' 
= -V P 2 v" 1 , 

n 2 3 


n 1 j 

j n +1 . 
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It is of course G^j which we wish to prove is exponentially 
damped. 

Let Gj be the Green's function associated with (5.4). 
We know 


(5.14) 


IG 1 ?! < K e ~ ak|n ‘ j| 

1 “ 


Now define = G ^ - G^j . If e V n corres_ 

ponding solution to (5. 3) , while x n is the solution 

to the inhomogeneous version of (5.4) , 


(5.15) 


x ,i - R x + kg„ 
n+1 n n 3 n 


Then u 
n 

(5.16) 


v - x is the unique bounded solution to 
n n 


u 


n+1 


R u + k 2 0(1) 
n n 


v 

n 


We thus have 


00 uu 

(5.17) u = k l H 1 * g x = k 2 l g" 0(1) 

11 j — -CO -J J j — -CO J J 


If we fix i and let - 0 f° r j ¥ i“l where 

z is an arbitrary vector in E^ , we then obtain, using (5.14) 
and the admissibility bound for (5.3), 


(5.18) IlHjII _< K 1 , for all n, i 

where is some constant independent of k. (The restric- 

tion "for small k" is always understood.) The use of (5.14) 
then implies the boundedness of G^ i.e. 

(5.19) "G*H < k 2 , 
where K 2 is independent of k. 
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We will now show that 


(5.20) UP^k) - P 1 II = O (k) . 

Observe that P 1 is independent of k. To obtain (5.20) we 

note that P^-P^ = ^o* ^ we ^ e ^ ne ^i-l = z, ~ 0 
(j ^ i) where z is an arbitrary vector we then have 
from (5.12) 

(5.21) v._ 1 = k G^T 1 g i _ 1 = k G?" 1 z . 

Substituting into (5.17) we obtain 


( 5 . 22 ) 


u. 


n 


kH?z = k 3 

1 



Gj 0(1) 



Now using (5.19) and (5.14) we see the sum in brackets in 

(5.22) is 0(l/k) whence we conclude that 

(5.23) IlH^t = O(k) 


and setting n = i = 0 we obtain (5.20). 

Note that this proves that the ranks of the stable and 
unstable manifolds are unchanged for small k. 

We can now show quite easily that (5.19) can be replaced 
by an exponential decay factor. 

For a certain small positive e consider new systems 
with the linear part of (5.3) and (5.4) multiplied by e i.e. 


(5.24) 


v 


n+1 


-ek fT , 
= e [R. 


n 


k 2 0(l)]vj 


(5.25) 


n+1 
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Now by Section 4 (see (4*9) f f - ) the inhomogeneous versions 
of (5.24) and (5.25) will both have (l£,l£) admissible, if 
e is small enough. But the new fundamental solutions are 


(5.26) 



-enk 

e 


V 


n 



(5.27) 



-enk 

e 


X 


n 



and clearly (5.25) will have an exponential dichotomy with 
projections P^ and P 2 if £ is small enough. 

Now if P^k) and P 2 (k) are the corresponding stable 
and unstable projections for (5.24) while Gj. s the 
Green's function for (5.24), the same analysis that led to 
(5.19) will yield for some constant K 3 , 

(5.28) llG j;e" - K 3 

In particular if j n+1, (5.2 8) reduces to., 


(5.29) 


II V 


n 


V. 1 ! < K_ 

1 — 3 


-ek ( j-nj 


and one part of the definition of an exponential dichotomy 
for (5.5) will have been shown if we can show P 1 = p 1 , 

t- 

P^ = p . But this is a simple consequence of (5.20). In 
fact if S+ , sf are the new stable. and unstable manifolds 
for (5.24) then (5.26) yields 


(a) s + C S + 

(5.30) 

(b) Sf C S_ 

But (5.20) shows that the ranks of and P^ are equal-' 
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(being equal to the rank of P-^ and similarly for P 2 and P 2 - 
This shows that the inclusions in (5.30) are equalities and 


that 


(a) 

(5.31) 

(b) 



Thus one part of the requirements of an exponential dichotomy 
for (5.5) has been shown and the other part will follow on 
replacing -e by +e. 

Lastly we point out that an exponential dichotomy for (5.5) 
implies one for (5.1) since in terms of fundamental solutions 
we have 

(5.32) W n = l" 1 (kA n )V n T(kA 0 ) 

whence the projections P^(k) , P 2 (k) for (5.2) are related to 
P and P 9 by a similarity transformation 

(a) P^(k) = T _1 (kA 0 )P 1 (k)T(kA 0 ) 

(5.33) 

(b) P^OO = T -1 (kA 0 )P 2 (k)T(kA Q ) . 
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6. Miscellaneous Results 


. This section is concerned with certain generalizations 
of the theory developed in the preceding section. We first 
remove certain restrictions which had been imposed previously 
in order to make the exposition clearer. In Section 6A we 
will remove the restriction that A be while in Section 6B 
we remove the restrictions (2,6) and (2.7) on the roots of 
the polynomial p(x) =0 which lie inside the unit circle. 

We then consider certain trivial extensions of the 
theory. In Section 6C we remark on the general (non-AP) 
case. In Section 6D we consider the convergence of the 
mean value of our sequence solution to the mean value of 
the solution of the ODE and finally in Section 6E we 
consider the weakly nonlinear case. 


- 49 - 



6A. Removal of the Differentiability Condition on A 


First consider the case when the homogeneous equation 


(6.1) y = Ay 
is exponentially stable; i.e. 

(6.2) H Y(t) Y -1 (s) II < K e " a(t ~ s > , t > s. 


where Y(t) is the fundamental solution to (6.1). We assume 

1 

the A is AP but is not C . Now Theorem 2 does not require 

2 

a perturbation 0 (k ) but is certainly valid for a perturbation 
k o(l) (k -► 0). This leads us to expect that the requirement 
that A be C 1 can be replaced by the uniform continuity of AP 
functions on the whole real axis . This is indeed the case 
as one can verify with little difficulty. 

In this subsection, we will merely trace through the 
proofs in Section 3 and indicate what changes must be made 
if A is not C 1 . 

Discretizing (6.1) we obtain the homogeneous difference 
equation 


(6.3) 


w . , = C w 
n+1 n n 


where C n is given in (2.11). Using the uniform continuity of 
A we see that (3.7) can be replaced by 


(6.4) w . = U (k A ) w + k o(l)w 

n+ 1 n n n 

where the term o(l) is uniform in n (this will not be stated 
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explicitly in the future) . If we also use the uniform 
continuity of A to obtain 


(6.5) T (k A - +1 ) = T(k A n ) + k o(l) 

we see that Theorem 2 will yield exponential stability for 
(6.3) provided we can show exponential stability for the 


n 


system 

r I+kA 0 

l6 ' 6) V n +1 = ( 0 " 

The lower block causes no difficulty and we are left with 
the system 


n 


(6.7) 


x 


n+1 


— [I + 


k A ] x . 
n n 


Finally we can show that (6.7) is exponentially stable by 

using Theorem 2 to compare the fundamental solution X n with 

Y (= Y (nk) ) . This requires some modification, as the proof 
n 

leading to (3.16) used the differentiability of A (see 


(3.13) ff). 

Following the notation of Section 3, we define 
(t,s) = Y(t)Y _1 (s). We can then write (compare with (3.14)) 


H 

(6. '8) H ( (n+1) k , nk) = I + kH(nk,nk) + k 

1 


d6 [H (nk+0k ,nk) 
- H (nk ,nk) ] 


= I + kA n + k | de [A(nk+0k) -A^] H(nk+6k,nk) 


+ k 


d0 A n [H(nk+0k,nk) - H(nk,nk)] 
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Finally we see that if we use (6.2) together with the 
uniform continuity of A on the whole real axis and the 
equation 

* 

(6.9) H(t,s) = A(t)H(t) 
then (6.8) can be rewritten as 

(6.10) Y n+1 = [I + kA n + k o(l)] Y n 

and Theorem 2 is immediately applicable. 

We thus have exponential stability for (6.3) and hence 

(A n ,A n ) admissibility for the inhomogeneous version of (6.3). 

Finally we point out that the results of Section 4, 

regarding admissibility in the case that (6.1) has a general 

1 

exponential dichotomy is equally valid if A is not C , as 
the contracting mapping principle (see (4.9) ff) would 
certainly be applicable if the perturbation is k o(l). The 
results in Section 5, however, use crucially the differenti- 
ability of A (in the argument involving the deduction of 
(5.18) from (5.17)) and it has not been possible to extend 
this result when A is not C^". 
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6B. Removal of Restrictions on the Roots of p(x) = 0 


In Section 2 we prescribed two conditions on the roots 
of p (x) = 0 which lie inside the unit circle; namely 

(6.11) x u distinct , u = 2,..., I, 

(6.12) x u ¥ ° / u = 2, .. 

(Recall that x-^ = 1 and | x^ | < 0 < 1 for u > 1) . We will 
now remove these restrictions. 

It should be clear by now that in the exponentially 
stable case (6.12) is unnecessary. In fact the fundamental 
solution to the homogeneous equation is used only for n >_ j , 
and the equation (3.5) shows that the invertibility of the 
linear term Cj_ is no.t required for n _> j . Of course the 
notation W^T 1 is no longer accurate, but except for this 
detail the results of Sections 2 and 3 are valid without 
assuming (6.12) . 

We deal next with the restriction (6.11). (Note that 
the root = 1 is always simple.) The removal of (6.11) 

in the finite dimensional case is trivial. In fact suppose 
the root x r has multiplicity q > 1; i.e. 

(6.13) p L (x r ) = ... = p q_1 (x r ) = 0 
p q (x r ) t* 0 . 

It is shown in Henrici (p. 214) that the effect of (6.13) 
is that the eigenvalue x r of U(0) (see (2.12)) will now have 
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nontrivial Jordan blocks. Thus U(0) can still be put in 
block diagonal form (see (2.17)) 

f 1 0 

( 6 . 14 ) U(0) = 

^ 0 B (0 ) J 

where B(0) is no longer diagonalizable, but has all of its 
eigenvalues bounded by 0 in the norm. Now since B(0) can 
be put in Jordan normal form with e instead of 1 on the 
superdiagonal, for any e > 0, we can certainly choose 
a matrix norm such that 

(6.15) I B (0 ) 5 < 6 


Since all norms are equivalent on a finite dimensional space, 
the proof of Theorem 1 can now proceed exactly as in Section 2. 

The case when E is infinite dimensional can be handled 
in exactly the same manner; however since we can no longer 
appeal to the theorems of linear algebra, we will have to 
carry out the proof in more detail. 

It is shown in Henrici (p. 214) , that in the scalar case 
a basis for the generalized eigenspace E r corresponding to 
the eigenvalue x is {e . } where 


(6.16) 



f x *' 1 } 
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fr U-r)x 2 - 'i ' 11 
r=l 

* • 

ff (n-r)x (n " {j_1)) 
r=l 

I » 

* 

ff (j - r) x 
r=l 

0 

0 


One can verify easily that in fact 


(a) U(0) e , = x e 1 

r , 1 r r , 1 


(6.17) 


(b) U (0) e . = x e . + 5 . . , 

r,j r r,] r,^-l 


j > 1« 


Also since the vectors e form a new basis the i x SL 

r 'J 

matrix with columns e . is nonsingular. Wow in order 

r,j 

to have e instead of 1 on the superdiagonal, we replace e r ^ 

by e . where 
1 r » 3 

(6.18) ■ e . = £-* 1 e . 

r I J ^ f J 

(see Bellman [9], p. 198). Here e > 0 is to be specified. 
The matrix with columns e^ ^ is still nonsingular but 
(6.17b) will now be replaced by 


(6.19) 


U(0) e 


r i J 


= x 


e . + 
r, 1 


e e 


r r j“l 
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Now the extension of this to infinite dimensional space 

is very simple and follows closely the procedure of Section 2. 

For any vector y e E define the vector e . (y) e E 0 by 

r r 3 


( 6 . 20 ) 




3 o-i ) 

TT U-rJx* 3 y 

r=l 

» 

j 

] T ( j+l-r) X y 
r=l 

0 


0 


Equation (6.19) clearly generalizes to 


( 6 . 21 ) 


U(0) e . (y) = x e . (y) + e e . , (y) , j > 

j-/j r r , j - 1 - / J 


u (°) e r ^(y) - x r e rfl (y) • 


We now define the space E . = (w e e o such that 

Z" 9 D X' 

w = e r j (y) for some y £ E>. Clearly the space E r ^ is 
closed for all r and j. Furthermore we have 


(6.22) 


E„ = 
l 


l E 


r t J 




This follows by the same argument as that given in Section 2 

r ?1 1 

(see (2.14) ff) . In fact if w = 


is a vector in E 


Z 


and we had an expansion 


(6.23) 


f X 1 } 

= T e . (y 

r,j J r,j 

x« J r '3 


V x A 


then the {y .} would be related to the {x. } by multiplication 
r , ] i 
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by a nonsingular i x l matrix and this establishes (6.22). 

It also follows from this that if we define a new norm on 

E. by II wfl u = max I y . It„ then Jwll u is equivalent to Hwll^ 
x, „ . r , 1 E 

r,j 

(max I x. II _) . 
i 1 

Now if E 1 is defined as 1 (the root x 1 is simple) 

while is defined as © l E . , then, exactly as 

r>l, j r,J 

in Section 2, we see that U(0) is in block diagonal form 
with respect to the decomposition E = E^ © E^, and we can 
write (see (2.17)) 


(6.24) 


U(0) 


■( 


D (0) 0 

0 B (0 ) 


where D(0) : -*■ E^ and is the identity, while B(0): E -*■ E' 

and for e sufficiently small we have from (6.19) 


(6.25) II B (0) II U < 9 

The proof of Theorem 1 can now be carried out exactly 

as in Section 2 if we make the final observation that the 

spectrum of U (0) is exactly (x u >. This should be obvious 

and follows from the fact that on each "generalized eigen- 

space" E = © l E r . , U(0) - XI (for \ ? (x u >) acts as 
r j r ' 3 

the matrix 

x r ~X e 0 0 

0 • 

(6.26) 0 

0 m . 

e 

0 . . . 0 x -X 

r 

and the inverse is obtained merely by inverting this matrix. 
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6C. Extension to the L Case 

— " ■■■ 1 oo 

Here we would like to make the simple observation that 
the difference equation theory, like the ODE theory, is 
essentially an theory which produces AP solutions when 
the coefficients are AP. 

In fact, if we now agree that derivatives are to be 

understood as instead of AP, the theory developed in 

Sections 3, 4 and 5 is entirely valid in the L ro case, with 

only one unimportant exception. The single exception is 

the proof given in (3.36) ff of the convergence of the 

f y n +£-i] 

solution w to the vector w_ = • in the case 

n ii [ • i 

l Yn J 

that the coefficients are not smooth. In fact if we 

refer to (3.36) and (3.37) we see that (3.37) need not 

be valid because the function y which is (we take A and 

f continuous) need not be uniformly continuous on the 

whole real axis . In this case we have not been able to show 

that w n -*■ w n uniformly for all n but only uniformly for nk 

lying in compact intervals. 

To see this let us write w„ (k) , w (k) to indicate 

n n 

explicitly the dependence on k. Now if the initial data 
for w R (k) is consistent, that is 

( Yl 1 

(6.27) w Q (k) - T ^-> I 

l Yl J 

for some vector y^ £ E, then it is well known (see 
Henrici, p. 244) that w n (k) must converge uniformly on 
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compact subsets, to some solution to the equation 

(6.28) y = A (t) y + f . 

But since every solution to (6.28) but one is unbounded, this 
limiting solution can only be the unique L ra solution y(t). Thus 
we must show (6.27) and the proof of this will in fact show 
that the vector y 1 in (6.27) is y(0). 

We assume at least that A is (with bounded derivative) 
so that the results of Section 5 are applicable. We can then 
write 

(6.29) w Q - w Q - k l g(k,(j~l)k) 

^ j — —CO -J 

where g is given in (3.37) , and satisfies the estimate 

-a k|n-j| 

(6.30) i! G“l! < K l e 

Observe that g(k,t) ■* 0 uniformly for t in compact intervals . Let 
T be unspecified for the moment and rewrite the sum in (6.29) as 

(6.31) w n -w n = k l G?g(k, ( j-l)k)+k l G°g(k, ( j-l)k) . 

|jk|<T 3 | jk | >T 3 

The first term -* 0 because g(k,t) -»■ 0 uniformly on the compact 

-a^T 

interval [~T,T]. The second term is 0(e ) as we can easily 

see from (6.30). Thus by first choosing T sufficiently large 
to make the second term small, and then choosing k small enough 
to make the first term small we see that IIWq-WqB -> 0 (k 0) and 
this establishes the convergence of w R to w n uniformly on 
compact intervals. 
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6D. Mean Value Properties 


Let a n e A n . 
the limit 


It is shown in Corduneanu (p. 48) that 


. n+N-1 

(6.32) iim - l a. = m(a ) 

N-*-» j=n 3 

exists and that this limit is uniform and independent of n, 
which we henceforth set equal to zero. This is of course 
the exact analogy to the mean value of an AP function 


T 

(6.33) lim ^ J f(t) dt 5 m(f) . 

T-voo 1 J 

Now if w n is our AP sequence solution, we would like to 

study the behavior of m(w ) as k -*■ 0 . Define w as 

n n 


(6.34) 


w n = 


y n-K-l 

* 

* 

1- Y n 


where y is the AP solution to the ODE 


(6.35) 


y = Ay + f 


Clearly 

(6.36) m(w n ) 


Now since (see (3.39)) 


m(y n ) ' 
m(y n ) 


(6.31) 


I w -w u 
n n °° 


o(l) 

k+1 


- 60 - 



we must compare m(y n ) with m(y) . We first claim for any 
AP function g we will have 

(6.38) m(g n ) ^ q- -> m(g) 

where g n = g(nk) . This is obvious because (6.38) holds if 
g is a trigonometric polynomial and for. arbitrary g, we 
simply approximate by a sequence of trigonometric 
polynomials. 

Now (6.38), together with (6.37) and (6.34), certainly 
yields 

iti(y) 
m(y) 

We now concern ourselves with the rates of convergence, 

p+ 

assuming smooth coefficients. If A and f are C (so that 
y is) we can replace (6.37) by (see (3.41)) 

(6.40) l * n “"n l « = 0(kP) 

where p is the order of accuracy of the scheme. It 
therefore follows that 

r m <Y n ) v 

(6.41) m(w n ) = • + 0(k P ) 

i «i(y n > . 

and it is only necessary to study the convergence of m(y n ) 

r+1 

to m(y) . In fact for any AP function y € c we have 
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(6.42) 


— T“+ 1 

m(y n ) « m(y) + 0(k ) 

Equation (6.42) is an almost periodic analogue of a 
theorem of Isaacson and Keller [10] for periodic functions 
(p. 340) and in fact the proof is very similar, and is a 
simple consequence of Taylor’s theorem. Letting t = nk 
and s € [0,1] we can write 

(6.43) y(t n +sk) = y n + (sk)y^ + ...+ - s ^, y^ + 0(k r+1 ) 


1 

where the remainder is uniform in t since y is AP. 
Integrating (6.43) from s - 0 to s = 1 we obtain 


'n+1 


(6.44) 


1 


n 


y(t) dt = k 


X 

I y(t n 


+ sk) ds 


= k tl ’n + l y n + ' • ’ (r+ITT y n + 0 'x" 1 ’ 1 


If we sum (6.44) from n = 0 to N , divide by Nk and let N -*• « 
we obtain 

— v 1 r r + 1 

(6.45) m(y) - m(y n ) + § m(y^) + ... ^jyy m(y£) + 0(k r+X ) 

Now if r = 0, i.e. y is the terms involving m Cy^ - ) 

through m(y n ) ar ® missing and we obtain (6.42). For r > 0 

we assume (6.42) holds for j = 0, ..., r-1. We can then 

. , 1 . _r 

write, as y is C , 

(6.46) 0 = in (y^) = m(y^) + 0(k r ) 
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r since 


and similar expressions for mfy 1 ) for i = 2,..., 
yi e c ( r- i) + l (we have used the obvious fact that the mean 
value of a derivative is zero), to obtain (6.42). If we now 
apply this to (6.41) we obtain 



6E. Weakly Nonlinear Equations 


Consider the equations 


(6.48) 

y = Ay 


(6.49) 

» 

y - Ay + 

f(t) 

(6.50) 

0 

y = Ay + 

h(t,y) 


If any pair of function spaces (B,D) is admissible for 
(6.49) the contracting mapping principle provides a technique 
to obtain D-solutions to (6.50) under certain conditions on 
the nonlinear term h(t,y), the most important being that 
h has a small Lipschitz constant (see Hartman [11], 

Chapter 12) . 

The case of (A, A) admissibility is particularly simple. 
We require that h be almost periodic in t, uniformly for y 
in compact subsets (see Hale [12] pp. 113 f f) . This simply 
insures that h(t,x(t)) will be AP for any AP function x. 

We also require that 

(6.51) llh(t,x^) - h(t,x 2 )ll oo _< e^lx^ - x 2 < 

for any x^,x 2 e E. The restriction on is simply 

(6.52) < 1 

where is the admissibility bound for (6.49) (see (0.3b)). 
The unique AP solution to (6.50) will simply be the limit 
of the iterates 
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(6.53) 


0 


0 

y = 

i+l 

y 

*i+l 

y 


is the unique AP solution to 
i+1 , . , . i, 

= Ay + h(t,y ) . 


This is of course an immediate consequence of the contracting 
mapping principle. 

The same proof will yield AP sequence solutions to the 
discretized version of (6.50). 


(6.54) w n+1 = C n w n + fc g(nk,k,w n+1 ,w n ) 


Here the nonlinear term is 


-1 


(6.55) g(t,k,w,w) = 


j- (a^I-k6^A(t+£k) [g^h (t,Xj_^) ^ 

l - 1 

+ l e h(t,x.) ] 

j=o J ■> 


( *A-l' 


where w = 


and w - 


*A-1 


. Observe that g is. 


for fixed k, AP in t uniformly for w and w in compact 
subsets, and that g will have an O(e^) Lipschitz constant 
which we call e, i.e. 


(6.56) II g(t,k,w 2 ,w 2 ) - g(t,k,w 1 ,w 1 )ll oo 

e max [lw' 2 - w^ll , Ilw 2 ~w^ll] . 

Now if K is the admissibility bound for the linear inhomogeneous 
version of (6.54) (see (3.30)) and y - eK < 1 then (6.54) 
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will have a unique AP solution , which in fact will be 

the limit of the iterates w , . defined by 

n / 1 


(6.57) w n;0 - 0 


w. 

w 


is the unique AP solution to 

C w . , , + kg (nk , k ,w , , . , w . ) 

n+l?i+l n n;i+l * n+l;i n;i 


This is an immediate consequence of the contracting mapping 
principle and we note for future reference that 


(6.58) 


l w -w 
n n; 


II 

m 00 


= 0(y m ) 


Now the first thing to consider is the convergence of 

f ^n+5,-1 


w to w defined as 
n n 


where y is the solution to 


(6.50) (assuming of course that (6.52) holds). In fact 

as was seen in Section 3 satisfies (6.54) up to an error 

k o(l) (0(k p+ ^) if everything is smooth). Hence the 
k^l 

difference z = w - w is the unique AP solution to 
n n n 


(6.59) z n+1 = C n z n + k[g(nk,k,w n+1 + * n+ l' w n + *n> 

- g(nk,k,w n+ ^/W n ) ] + k o(l) (0(k^ +1 )) 

and we have 

(6.60) II z II < yll z II + o(l) (0(k p )) 

n 00 — n 00 

and since y < 1 this settles the question of convergence. 

We now examine the question of the computability of 
the solution in the exponentially stable case? that is when 
we have constants K and a such that 
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(6.61) II W wT 1 !! < K e “ ak < n “^ n > j, 

n 3 — ~ 

where W is the fundamental solution to the homogeneous . 
n 

version of (6.54). 

To do this we use (6.58) to reduce the problem to a 

study of the effect of round-off errors and errors in initial 

data on the solution w for a fixed m. This will be a 

n;m 

simple consequence of the stability results of Section 3 

(see (3.42) ff) for the inhomogeneous equation. 

In fact the results of Section 3 show that at the 

first iteration we would solve (for n 0) for a solution 

w , which is related to w„ , by 

n ; 1 n ; i 

(6.62) w* . = w , + 0(k r ) + O (e akn ) . 

n ; i n j x 


We now consider the effect of this error on the second 

iterate. Assuming/ for the moment, exact calculations and 

* . 

exact initial data, we obtain a solution w n -2 which solves 
the equation 


(6.83) w n+1 . 2 - C n w n;2 + k g (nk,k,w n+1(1 ,w n!l > . 

If we let z n;1 - w* ;i 


- w for i = 1,2 then we obtain 
n ? i 


< 6 - 64) z n+l;2 - C n z n;2 + k [ » (nk ' k ' w ntl! 1 + Z n+1; l' w n ; l +Z n;l ) 

- g(nk,k,w n+1;1 ,w n;1 )] 


The last term can be written (using (6.62)) as 
k 0(k r ) + k 0(e~ akn ) and we have 
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(6.65) z „ = k l W wT 1 0(k r ) + k l W W" 1 0(e" akj ) 
n,z j=l n 1 j=l n 1 


r -akn r _a i kn 

= 0 (k r ) + 0{e a kn) = 0(k ) + 0(e x ) 


where can be arbitrarily close to a. We have used 

(6.62), (6.61) and the trivial fact that t = 0(e < ^ t ) for 


any positive 6. 


It is now obvious that (6.65) is preserved if we 

* 

include in the equation for w ~ (see (6.63)) the effect 

n ; z 

* 

of roundoff and errors in the initial data of w^ . 0 

n ; z 

(compared with the initial data for w -) and it is also 


obvious that this process can be repeated for m itera- 


tions , for any fixed finite m, and this establishes the 


computability of the solution w^. 
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Part II. Extensions to Partial Differential Equations 


We will now use the techniques of Massera and Schaffer/ 
in the simple exponentially stable case/ to obtain AP solu- 
tions to a certain class of inhomogeneous partial differential 
equations (PDE's) to which one can apply ODE formalism by 
use of the Hi lie Yoshida Theorem. The homogeneous equation 
has been extensively studied by Krein [1 3] , to which we will 
refer often. A more succinct study of the Hille-Yoshida 
Theorem can be found in the appendix of Lax and Phillips [14J . 

Throughout Part II E will denote an infinite dimensional 
Banach space. 
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7. Existence Theorems 


Consider the linear equation 

(7.1) y = By 

where B is an unbounded operator. The study of the solutions 
to (7.1) has been answered by the Hille-Yoshida Theorem in 
the case that the resolvent R^(B) (- (AI-B) "*”) satisfies 

the inequality 

(7.2) II (B) II <_ j , A > 0 . 

A more general formulation is given in Krein (in particular 


Section 2, 

Chapter 1) . Krein shows that if: 

(7.3) 

(a) 

B is closed 


(b) 

B is densely defined 


(c) 

ffR^(B)B < , R.P. A > to , 

A (Re A- to) n 

then 

there 

exists a semigroup , t > 0 , such that: 

(7.4) 

(a) 

Z t is strongly continuous 


(b) 

ESI 

O 

II 

H 


(c) 

II Z t l! £ Me wt 


(d) 

lim [Z^-Ije/A converges iff e e D g 
A-*0 

and in that case it converges to Be 


(e) 

if e e D 0 then for t >_ t^ the unique 


solution to (7.1) such that yft^) - e 

is simply y(t) = Z e (uniqueness is 

t-t 0 
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shown in Theorem 2.7, page 47) 

(f) z. : D„ D- and commutes with B for t >_ 0 . 
t B B 


For a well posed problem we obviously require w <_ 0 in (7.4c) . 
To insure exponential stability we are going to proceed as 
if a) as o and introduce our own damping term which will be 
allowed to depend on t. The reader should note that if 
w < 0 this term is not necessary. It is well known (see • 
Krein, page 43) that one can always introduce an equivalent 
norm so that the factor M in (7.4c) can be replaced by 1, 
however we will eventually approximate B by a family of 
bounded operators where 

h B t 

(7.5) I' z£ = e n I £ M 

and we may not be able to introduce a norm so that M can 
be taken as 1 for all h, and thus we will leave M unspecified. 

Now to make the homogeneous equation exponentially 
stable we introduce a damping term -6(t) and consider 
the equation 

(7.6) y = [B - <5 (t) ]y 

For conditions on 6 we first let it be a scalar AP function 
such that 


(7. 7) 


R.P. (m(<5) ) ' > 0 


This implies that if g(t,s) = exp-j- dr 6(r) j- then 


h: 


(7.8) 


|g(t,s)| £ K e , S ^ , 


t > s , 
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as can easily be seen. 

It is also possible to permit 6 to be an AP operator 
function satisfying certain commutivity properties. Speci- 
fically if Y (t) is the fundamental solution to 

{7.9) Y = - oY 

then 

(7.10) (a) 8 Y ( t) Y _1 (s ) II < K e _a(t “ s) , t >_ s 

(b) 6(t), Y (t) , Y ^(t) commute with 

for all t and all r > 0 

(c) 6 (t) , Y { t) , Y ^(t) map D d into itself 
and commute with B. 

These hypotheses permit us to essentially treat 6 as a 
scalar in the following analysis. If we define 

(7.11) g(t,s) = Y(t)Y~ 1 (s) 

then the "fundamental solution" to (7.6) is 

(7. 12) g (t,s) Z t _ s , t > s, 

and 

(7.13) Ug(t,s)z I < KM e -0t(t “ s) = K e ~ a ( t-s ) 

u"S 

after a redefinition of K, 

Now consider the inhomogene ous equation 

(7.14) y = [B-6 ]y + f(t) 

where f e A(E) or C m (E) . We want to study bounded (for all t) 
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solutions to (7. 14) . First we see that any such solution 
is unique. In fact if w(t) were a bounded solution to 
(7.6) we would have (for t >_ s) 

(7.15) w ( t ) = Z t _ s g( t,s) w(s) 

and w = 0 follows immediately on letting s -*■ 

It is equally simple to see that any solution must 
be given by 

t 

(7.16) y (t) = | g(t,s ) Z t _ s f(s) ds = C(f) . 

—00 

This is in fact an immediate consequence of letting t Q -*■ - 00 

in the variation of constants formula 

t 

(7.17) y (t) = g(t,t Q )Z t _ t y(t Q ) + | ds g (t/S) Z t _ g f (s) 

t 0 

which is proved in Krein (Theorem 6.1/ page 129) for the 
autonomous case, and it is a trivial matter to see that the 
same proof will work if 5 is time dependent. Note that (7*17) 
need not be a solution to (7.14), but any such solution 
must be given by (7.17). 

We now consider the operator C defined by (7.16). 

C is a bounded operator mapping L ffi L w . Boundedness is 
an immediate consequence of (7.13) and in fact 

(7.18) iy'„ 5 k (- 

The integrand is obviously continuous if f ^ and the 
only point in question is its measurability for a general 
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f 6 L # (for measurability of functions in a Banach space 

see Hille and Phillips [15]). However it is shown by Krein, 

and will be used crucially very shortly, that Z is the 

B r r 

strong limit, pointwise, of e n where B is a bounded 

n B r 

operator (Theorem 2.9, page 48). Since e n will be 
norm continuous this establishes the integrand as the 
pointwise limit of measurable functions, hence measurable. 

We will now show that y as defined in (7.16) is a 
solution to (7.14) . We will make use of the above mentioned 
approximation theorem of Krein (Theorem 2.9). It is shown 
that there exists a squence of bounded operators B n such that 

(7.19) (a) I (B-B n )el ~ n ~ > 0 (Ve e D fi ) 

B r 

(b) II (Z -Z^)ell - ■ > 0 Ve , r > 0, where Z^, = e n 

jl t n~ K3 ° — t 

(c) II Z^l _< M , r :> 0 

(d) B^, Z^ (r >_ 0) satisfy the same commutivity 

relations with <$(t) (7.10b,c) as does B and Z^. 

We only have to point out that (7.19d) follows from the 

2 

explicit definition of B as -A I - A R. (B ) where X -*■ 00 

n n n A n n 

(see Krein, page 49). We point out that the operators B^ 
cannot in general be considered as spatial discretizations. 

A family of such approximations will be assumed in the 
following section. 

Now consider the operators C n (f) defined by 
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t 

(7.20) y n - j g(t,s) z£_ s f(s) ds = C n (f) . 

— CO 

It follows immediately that the operators are uniformly 
bounded operators mapping L m into and in fact 

(7.21) iy n ' B i Kiifi'cc 

where we take the same bound as in (7.18) . It is also 
immediate that y^ is the unique solution to 

(7.22) y n = [B n -S]y n + f 
and y is AP if f is. 

We are going to show that C strongly in a certain 

subspace of L^. Specifically define S to be the subspace 
spanned by the functions whose range has compact closure. It 
is well known (see Amerio [16]) that ACS. For f £ S we 
will have 

(7.23) C n (f) * C(f) 

To show (7.23) note that by the. uniform boundedness of C n 
it is sufficient to show it for a dense subset of S and by. 
the definition of S it is sufficient to take f as 

n 

(7.24) f = r 

j=l 3 3 

where e. e E and Xj is the characteristic function of a 
measurable set. By linearity we only have to work with 
one term in the sum in (7.24) . Thus if I is some measurable 


- 75 - 



set let f = x x e • Then w (t) " y(t) - Y n (t) is explicitly 

OO 

(7.25) w n (t) = | x( t “ r ) g(t,t-r) [Z r ~Z^]e dr . 

0 

Estimating (7.25) we obtain 

00 

(7.26) [lw n“co 1 K I e " ar,[ (Z r -Z”)e# dr 

0 

and (7.23) is an immediate consequence of the Lebesgue 
dominated convergence theorem. 

Observe that (7.23) implies that C(f) is continuous 
if f 6 S and it is AP if f is. 

Now consider the almost periodic case. First we point 
out that if <5 and f are C 1 (AP derivatives) then y is 
also and in fact 

(7.27) y = C(-6y + f) 

This is trivial if we observe that 

(7.28) y n = C n (-6y n + f) 

and if y is defined to be the right-hand side of (7.27) we 
have 

(7.29) y-y n = !C-C n ](-6y) + C n (-5 (y-y n ) ) + Ec-c n l(f) , 
and clearly 

(7.30) 8 (y-y n )i TO 0 

Equation (7.27) is now an immediate consequence of letting 
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n °o in the equation 


(7.31) y n (t 1 ) - y n (t 2 ) = | y n (t) dt 

We next show that y satisfies an integrated version 
of (7.14). Specifically 


(7.32) y(t x ) - y(t 2 ) 


t 2 Z 2 t 2 

= B | dt y(t) - | dt S(t)y(t ) + j dt f(t) . 


Since Dr, is dense and B is closed it is sufficient to show 

(7.32) when f is of the form 


(7.33) f = e lXt b , b S d B* 

Now observe that Bf is AP. It certainly follows that 


(7. 34) By = C (Bf ) . 

In fact this holds whenever Bf e and is a simple conse- 
quence of the closure of B, approximating the integral in 
(7.16) by Riemann sums over finite intervals and then 
noting that B applied to each sum is the Riemann sum 
approximating C(Bf). 

Now the same property certainly holds for B n and C n i.e. 
(7.35) Vn = c n (B n f) ‘ 

Finally since by its construction B n commutes with 
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(Krein, Theorem 2.9) we can write 

(7. 36) By = C (B f) . 

n n 

Since y satisfies (7*22) we have 
u n 


(7.37) y n ( tl ) - y n (t 2 > 


- B n f dt y n (t) * 


dt 6{t)y n (t)+ f dt f(t) 


U 

/ 


= B dt y ( t) + [B 


‘n' B1 f 


dt y (t) + B 


n 


d t (y n ~y) 


2 ^2 
- | dt 6{t)y n (t) + J dt f(t) . 


1 ''l 

2 


where the fact that j ^ dt y(t) e D g is a consequence of 

^1 

(7.34) and the same type of closure argument as followed 
that equation. 


If we now let n •* <® in the last equation in (7.37) we 

ft 2 

dt (y -y) - > 0 because 

fc l 


see that the term B 


n 




(7.38) B n (y n -y) - t c n ~ Cl (B n f) = [C n -C](Bf) + (C^C] ( (B^B) f ) 


and the right-hand side of (7.38) clearly + 0 in the sup norm 
as n ■+ 00 . It therefore follows, on letting n -*■ 00 in (7.37), 
that (7.32) does in fact hold for functions of the form (7.33) 
and thus for all f e a. 

Now if 6 and f are (hence y is) we can show very easily 
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that (7.32) implies (7.14). In fact since y is C 1 we can 

__ t+h t+h 

t+h t f 

B I ds y (s) J ds 6 (s)y (s) J 


let h 0 in the equation 

t+h 


f (s) 


(7.39) 


y (t+h) -y (t) 


where h can be positive or negative, using the closure of B. 
Similarly (7.14) will hold if we only require Bf is AP. 

In fact since in that case By is AP, B can be brought inside 
the integral in (7.32) , which can then be differentiated 
directly. 

We point out that in both these cases both y and By are 

also AP. For future use we point out, using the integral 

2 

equation (7.32) applied to (7.27), that if 8 and f are C 
2 

then y is C and y satisfies 

• * 

(7.40) y = [B-6 3y - 6y + f 

i.e. (7.14) can be differentiated formally. (Note in parti- 
cular that if Bf is AP, B 2 y is AP. ) This argument can 
obviously be extended; namely if f and 6 are C P then y is C p 
and we have 

(7.41) y r = [B-6]y r 1 + f r 1 - (Sy) r ^ + <$y r ^ r=l,...,p. 

We will now examine briefly the case that f e s. 

• • 

We still have y n -*■ y but (7.30) need not hold since Sy and f 
need not belong to S. However we see very easily that if 
f e L ra and 
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x = C (g) 

(7.42) 

X n ^ C n ( 9> 

then we will have x n (t) ♦ x(t) pointwise and boundedly. In 

fact if w n (t) ® x(t) - x n (t) , then 

00 

(7.43) llw n (t)S <K | e“ ar n (Z r -z") g(t-r)ll dr 

0 

and the Lebesgue dominated convergence theorem is still 
applicable and hence (7.31) will still yield (7.27), if we 
use the bounded convergence theorem and the fact that 

(7.44) y(t)-y (t) = [C-C ] (-6y+f ) (t) + C (-6(y-y ))(t) > 0 

ii n n n n-^ 00 

Thus y has a bounded derivative. Now if we assume that f is 
of the form 

(7.45) f = Xj b b e d b , 

then the derivation of (7.32) from (7.37) and (7.38) still 
follows, if we use the bounded convergence theorem for inte- 
grals over finite intervals. Since finite sums of such f's 
are dense in S, (7.32) holds for all f e s and if Bf e S or 
6 and f are C^, y will still be the unique C m solution to 
(7.14). Thus the theory will hold for a certain class of 
bounded, non AP, inhomogeneous terms, but is not as simple 
as the AP case and in the future we will restrict ourselves 
to this case. 

Finally we point out that (7.32) can be derived by 
brute force dif ferentation of the formula (7.16) but use of 
the approximating operators B n permits a cleaner and more 
straightforward development. 
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8. Discretization in Space 


The operators B n which proved so useful in the 
preceding section need not correspond to a discretization 
in space. We therefore postulate the existence of a 
family of bounded operators , where h will generally 

play the role of a spatial grid size, such that: 


(8.1) (a) 

(b) 

(c) 


(d) 

(e) 


(f) 


B, are defined and bounded for h £ (O/h^l 
h u 

I B. II = 0(-“) , m > 0 

h h ra 

B^ generate bounded semi- groups (i.e. 

B, t h 

H e h = Z^ll £ M for t £ 0) 
e e D fi implies II (B-B h )ell h _^ -> 0 

If <5 is not a scalar function, then the 

commutivity relations (7.1Q)b,c hold with 

h 

B and Z r replaced by B h and Z r 

If e £ D d (or any dense manifold in D 

invariant under Z^ for r >_ 0) we have 

lB,_Z ell < K (I ,e) for r e i, where I is any 
hr — 

compact interval on the nonnegative real axis. 


We point out that (8. If) is designed to insure that 

Z h -> Z strongly, for all r > 0. This is certainly reasonable 
r r ~ 

since B^Z^e converges to BZ r e. It is of course trivial if 
B commutes with Z . Also note that the convergence in (8. Id) 

h r 

does not deal with the order of the approximation of B h to B. 
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This will be dealt with by a subsequent assumption. 

Before proceeding further, we indicate a simple example 
of the case we are considering. 

(8.2) (a) E = L 2 (0,2 tt) 

£ 

(b) B = -r— (with periodic boundary conditions) 

, u(x+h) -u(x-h) 

(c) V 2h 

In this example (8. If) is unnecessary as B^ in fact commutes 
with Z^. 

Now associated with B^ we have the equations 

(8.3) y = [B h -6]y 

(8.4) y = [B h -6]y + f 

The homogeneous equation (8.4) has as its fundamental solution 

B.t 

Y^(t) = g(t,s) e , and is exponentially stable, i.e. there 
exist constants K and a independent of h such that 

(8.5) II Y h ( t ) 1 ( s ) II < K e” a (t “ s) , t > s. 

Associated with (8.4) we form the operator where 

t 

(8.6) y h = C h (f) = | g(t,s) z£_ s f(s) ds 

— oo 

for any AP f, y^ is the unique AP solution to (8.4). 

Equation (8.5) shows that the operators are bounded 

uniformly in h, i.e. 

<8-7) 1 y h > „ i Kl £! „ 

where K is independent of h. 
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Now our first task is to show that C strongly. 

Using the results of Section 7 (see (7.23) ff) it is suffi- 
cient to show that z“ z strongly, for r >_ 0. This will 
be an easy consequence of (8. If). In fact if e e Dg (or 
the manifold described in (8. If)) and if y = Z^e , 
y^ = Z^e , then if w* 1 = y-y^ we will have 


( 8 . 8 ) 


« h = B h w h + t B - B h ] y 

w h (0) = 0 . 


From (8.8) we can write 

t 

(8.9) w h (t) = | ds z£_ s [B- V Z s e 

0 

and the fact that Bw h (t)i 0, in fact uniformly on compact 
intervals, is an immediate consequence of (8. If) and the 
Lebesgue bounded convergence theorem. The result for any 
e e E follows from the uniform boundedness of the 
operators Z* and the denseness of the manifold of (8. If). 
Now the fact that 


(8.10) 1 (Y-Y^co “hTo _> 0 

for all f e A (in fact for all f e S) is not any improvement 
of the corresponding result for the operators of Section 7. 
We want to show that the difference in (8.10) approaches zero 
as fast as some power of h, and to do this we must require 
that B^ approximate B up to a certain degree of accuracy. 
Specifically we suppose that there exists an operator L 
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(unbounded) with D L £ D g and such that if e £ D L we will 
have 

(8.11) 11 (B-B^) ell £ C h ^ II Le H where j > 0 

and C and j are independent of h and e. 

If we have Lf(t) AP we would like to show that (8.11) 
implies 

( 8 . 12 ) i (y-y h )i„ = 0 {h j ) . 

Unfortunately (8.12) requires a further assumption; namely 

(a) Z r : D l - D l 

(8.13) 

(b) If Lf is AP then llLZ r f(t)H will be 
uniformly bounded for r >_ 0 and all t. 

Assumption (8. 13) is of course trivial if L commutes with Z r « 
It is now a simple matter to show (8.12). In fact 
w^ = y-y^ is the unique AP solution to 

(8.14) w h = [B h -S]w h + (B-B h )y 
whence by (8.7) we have 

(8.15) 11 w h 11 oo = °(“ [B-B h ]yllJ 

and so it is only necessary to study (B-B^)y. But now if 
Bf £ A, which will generally follow if Lf £ A since L will 
be a higher order operator than B, we can write 
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(8.16) [B-B h ]y - J dr g ( t, t-r) [B-B h ] Z r f (t-r) 

0 


and clearly 


(8.17)' II {B-B h ]yll oo = 0(h j ) 

is an immediate consequence of (8.13). 


C\ 
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9 . Discretization in Time bv the Multistep Linear Method 


We will now apply the theory developed in Section 3 
to the spatially discretized equation 

(9.1) y h = [B h -<S] + f • 


Since is bounded the theory is immediately applicable and 
we can assert that for k sufficiently small {but depending 
on h) there will be an AP sequence w n (k f h) which will converge 
to w n (k,h) defined by 


(9.2) 


w n (k ,h) 


r y h ( (n+&-l)k) ' 

m 

y h (nk) 


and if everything is smooth in t we will have 


(9.3) 


I W n - W]1 l 


‘ °h (k£, » 


where the bound in (9.3) may depend on h. It is only left 
to show that the bound in (9.3) is independent of h and to 
discuss the relationship of k to h. We will show that what 
is required is 


(9.4) 


k + 


k 

h 2ra 


- 0 ( 1 ) 


i.e. there exists a constant c independent of k and h, 

2 m 

such that if k + k/h < c the difference equation corres- 
ponding to (9.1) will be exponentially stable with constants 
independent of k and h. 
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The quadratic dependence on h in (9.4) is restrictive 
and will be removed in the following section. 

Throughout this section we will use the notation of 
Sections 2 and 3. We will also take 6 to be C^. 

We now fix h and consider the homogeneous difference 
equation on (see (3.7)), 

(9.5) w n+1 = U(kB h ~ k6 n ) + k 2 0(l)w n . 

It is clear that this will hold uniformly in k and h if 
k + k/h m is sufficiently small, i.e. if we are in a region 

(9.6) k + = 0(1) . 

h 

Of course restrictions of the form (9.6) are included in. (9.4) 
Also note from the explicit definition of the full linear 
term in (9.5) (see (2.11)) the last term in (9.5) will be 
uniform in k and h. For the rest of this section this will 
be understood for all "O" signs unless stated otherwise. 

If W n (k,h) is the fundamental solution to (9.5) it is 
clear that all we must show is that in a region of the form 
(9.4) we will have 

-a k(n-j) 

(9.7) BW n W j * - K 1 e ' n - 3 ' 

where the constants K 1 and are independent of k and h. 
Following the procedure of Section 3 we let 

(9.8) v „ " T(kE h ■ kS h )w n 
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and we then have in a region of the form (9.6) 

(9.9) v n+1 = L(kB n - k« n )v n + k 2 OU)v n 

Using the structure of the block diagonal operator L(z) 
derived in Theorem 1 we can write this as 


(9.10) 


v 


n+1 


[I+k(B -6) 
n n 

0 


0 

B 


n 


v + kO(k + -£-)v 
n li m n 


Here we have used the notation B n for the lower block (see 

(2.26)) so as not to confuse it with the operator B^. 

Observe that in (9.10) we have the first appearance of 
2m 

the term k/h . This is because in order to apply Theorem 2 
we must separate out a factor of k in the perturbing term. 

Now it certainly follows that in a region of the type 
(9.6) we have 

n-1 

(9.11) I TT B. II < K 6 n 3 , n > j, 

i=j 


where and 0 are independent of h and k. We therefore 

only have to consider the upper block in the leading term 

of (9.10) and in order to apply Theorem 2 we are going to 

compare this with Y, v = Y, (nk) where Y, is the fundamental 

n ;h n h 

solution to 


Y h - iv (|I h 


(9.12) 

Observe that the estimate 

(9.13) If y .y“* ,11 < K e " ak ( n -j) 

n;h ;j;h - ' 


n > j 
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holds, where K and a are independent of h. 

We will show that Y , satisfies the difference 

n ? h 

equation 


< 9 - 14 > W -■ t I+kD n tk - h)IY n;h 

= [I + k(B h -« n )]Y n . h + kO(k + ^;)Y n . h • 

Observe that the last term in (9.14) is of the same order 
as the perturbation in (9.10). 

Equation (9.14) follows easily from the same argument 
as in Section 3 (see (3.13) and (3.14)) if we simply observe 
that 


(9.15) 


(B h _<S) = “ 6 


and that (3.15) holds in this case for t > s by (8.5). 

We can -now apply Theorem 2 comparing the equation 

(9 - 16 > 5 n+l;h * II+k, W ,5 n,h 

with (9.14). We simply note that the proof of Theorem 2 
indicates that the terms e Q , a ^ and K 1 (see the statement 
of Theorem 2) depend only on the constants K and a in (9.13) 
and are thus independent of both k and h. Another applica- 
tion of Theorem 2 comparing the equation 


(9.17) 


V 


n-f 1 


■( 


I + k(B h -6 n ) 


V 


B 


n 


n 

with (9.10), together with the boundedness of the operators 
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T(k (Bj^-6 ) ) , T ^ (k (B h -6 n ) ) -*- n the re 9i° n described by (9.6) 

2 ^ 

will yield (9.7) with the restriction k + k/h sufficiently 
small (i.e. (9.4)). 

If we now refer to the proof of admissibility in the 
stable case given by (3.31) ff we see immediately that (9.7) 
implies (L^,L^) admissibility for the inhomogeneous version 
(9.5), with an admissibility bound independent of k and h. 

It is only necessary to show that the convergence of w R to w n 
is independent of h (see (9.3)). This will follow from (3.40) 
if we can show ly^ + ^8- oo will be bounded uniformly in h, 
assuming that <S and f are smooth. 

In fact if 6 and f have p+1 derivatives the differenti- 
ability properties of y^ are the same as those of y n introduced 
in Section 7. In addition, the formula given in (7.28) can 
obviously be extended (for the operators C^) to yield 

(9.18) y P+1 = C h (f P+1 + (-<$y h ) P+1 + <Sy£ +1 ) 

(note that the argument of in (9.18) does not involve y P+ ^) 

and the uniform boundedness of the operators immediately 

yields the uniform boundedness of the derivatives of y, . 

_ y ( (n+£-l) k 

If we now define w = : 

{ y (nk) j 

(uniformly in h) together with (8.12) implies 

(.9.19) Iw -w II - 0(k p +h j ) = 0 (h 2mp +h^ ) 

xi n 

where we have used (9.4) for small k. 


, then (9.3) 
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10. The Lax-Wendroff Scheme 


2 m 

The term k/h in (9. 4) can be improved for operators 
such that the Lax-Wendroff scheme is stable for the homogeneous 
equation 

(10 . 1) w = Bw 

We will then be able to obtain AP sequence solutions under 
the more favorable restriction 

(10.21 k + £- = OC1) . 

h m 

First of all we require that 6 and f are C 2 . As shown 
in Section 7 this implies that y is C 2 « It implies further 
that y. satisfies 

(10.3) y = !B-6]y - Sy + f 
i.e. that the basic equation 

(10.4) y = [B-6 ]y + f 

can be differentiated formally. 

We can now apply the Lax-Wendroff scheme to (10.4). 
Proceeding formally we write, assuming Bf is AP, 

k 2 .. 

(10.51 ylt+kl = y(t) + ky(t) + — y (t) 

k 2 2 

= y(t) + k[B-6 ]y ( t) + IB-6] y(t) 

k 2 * 

+ j- (~<Sy + (B-6 ) f +fj + kf ( t) = 
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= [I+kB + y ~ B 2 ]y (t) - kfiy 
+ k[-|- 6 2 - | 6 - k6B]y+ k[f+ |(B-S)f+ | f] 

We have written this as an inhomogeneous difference equation 

in order to apply techniques already developed. 

Now in the expression for y we approximate B by B^ 

2 2 

and in the expression for y approximate B by B ^ , for a 

certain constant c > 0 to be determined (c = 1/2 for the 

example described in (8.2)). Of course we will have to 
2 2 

assume that B^ approximates B in the same sense that B^ 
approximates B, but we will leave the precise hypothesis 
for later. 

Applying these approximations to (10.5) and letting w n (k,h) 
stand for the dependent variable we obtain the difference 
equation 

(10.6) w n+1 (k,h) = U(k,h)w n ~ k(S n w n + k R(n,k,h)w n + kf n 
where 

k 2 2 

(10.7) (a) U (k , h ) = I + kB h + ~ B c h 

< b > V f n + I [B ch ‘ 6 n lf n + I *n 

(c) R(n,k,h) = I - I «n ' kS n B ct> = 0(k+ ^ ' 

Now U(k,h) represents an approximation to (10.1). We 
choose c so that for 

( 10 . 8 ) = 0 ( 1 ) 
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this discretization is stable, i.e. 

(10.9) B U(k,h) r ll < H , r >_ 0. 

We are going to use Theorem 2 to show that if UR(n,k,h)li 
(i.e. k + k/h m ) is sufficiently small, then the homogeneous 
version of (10.6) is exponentially stable; that is if 
W n (k,h) is the fundamental solution to the homogeneous 
equation, then there exists constants and ot-^ independent 

of k and h such that 

-a.k (n- j) 

(10.10) HW^W^II < K x e , n > j, 

where for simplicity we have suppressed the dependence of W r 

on k and h. Once we have obtained (10.10) , (L^,!^) (and 

(A ,A ) ) admissibility will follow in the usual manner, 
n n 

To obtain (10.10) we observe that since we can allow 
k + k/h m to be as small as required. Theorem 2 permits us 
to neglect the perturbing term R. Also by neglecting a 
pQ;E"turb ation of the same order as kR we can work with the 
simpler equation 

(10.11) w^ +1 = [I-k6 n ] U(kh)w* . 

Now consider the function g(t,s) defined in (7.11). 
Using the same argument given in Section 3 (see (3.16)), we 
see that g n = g(nk,0) will satisfy the equation 

(10.12) g n+1 = (I+kD n lg n = [I-k« n +k 2 o (1) )g n 

where the "0" sign in (10.12) is independent of n. Thus by 
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again neglecting a perturbation of the form kO(k+k/h m ) 
we obtain finally the equation 

(10.13) w n+1 = [I+kD n ] U<k,h)w n . 

But the fundamental solution to (10.13) is simply 

_-l 

(10.14) W n Wj = g(nk, jk)U(k,h) n 3 

as one can verify immediately, using the commutivity 
properties of g n with (see (8.1e)) in the case that 

6 is an operator. Now the estimate 

-1 _ -ak (n- j ) 

(10.15) ! WW j I < K e n > j, 

follows from (7.10a) and (10.9) and we can now apply 
Theorem 2, noting, again that the terms K^, a^, £q depend 
only on the terms K and a (see the hypothesis of Theorem 2) , 
to obtain (10.10). 

We now have (L ,L ) (or (A ,A ) ) admissibility for 
n n n n 

(10.6) and it is only necessary to study the convergence in 

the case that the inhomogeneous term is given by (10.7b). 

2 2 

To do this we must now assume that approximates B ; namely 

that there exists an operator L such that if e e n d 9 

L 

then 

(10.16) II (B 2 -B^)ell < C 2 h 3 llLell 

We also assume that L satisfies the same assumption as L 
(see (8.13)) so that if Lf and Lf are AP we can conclude 
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(using the same argument as in (8.16) and (8.17)) 


(10.17) (a) I (B-B h )yH ^ = Ofh 1 ’) 

(b) II (B 2 -B^)yll oo = 0(h j ) . 

If w n (k,h) is the unique AP solution to (10.6), we can now 

show that converges to y R uniformly inn. 

3 

In fact since y is C , y n satisfies (10.5) up to an error 
0(k 2 ). If we then make use of (10.17) we see that y n will 
satisfy (10.6) with an error kO (h^+k ), where this is uniform 
in n. Then as the admissibility bound of (10.6) is independent 
of h (by (10.10)), we can conclude 

(10.18) llw^y^ = 0< h:j+k2 ) = 0(h j +h 2m ) 

where we have used (10.2). 

For the example described in (8,2) we will have 

g 3 g 4 . 

(10.19) L - — j and L = , c * j , m = 1, j = 2, 

3x 3x 

(periodic boundary conditions) 

Finally we would like to point out that in the term f n 
in (10.7b) one can replace B ch f n by Bf, without changing 
the order of the error. This is simply because Lf is AP, 
hence 

(10.20) II (B-B h )fl w = 0(h j ) . 
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11. Linear Perturbations 


We now consider the perturbed system 

(11.1) y = [B— <S ]y + D ( t) y + f(t) 

Here D(t) is a bounded operator for each t and is almost 
periodic as an operator function. We also require 

(11.2) BD^ £ Eg 

where Eg is to be specified. 

The first step in finding AP solutions to (11.1) is 
to solve the equation 

(11.3) y = C (Dy) + C(f) . 

This equation can be solved uniquely by the contracting map- 
ping principle provided p. = EgK < 1, where K is the norm 
of C (see (7.18)). The contracting mapping principle also 
shows that the assignment of y from f is a bounded operator 
and in fact 

U1.4) lyl„ £ = Kjf'^ 

where K-^ depends only on Eg and is independent of D. 

We want to show that y satisfies (11.1), and if we 
assume that D, 6 and f are c\ the results of Section 7 
show that it is sufficient to show that y is C"*". This will 
follow if we observe that y will be the limit of iterates 
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0 


y^ defined as 

(a) y 0 = 

(11.5) 

(b) y i+l = c(Dy i ) + C(f) 
We will then have for each i 


(11.6) y i+1 = CtD^) + C (-6y i + Dy.^ + f) 


If we then define y as the unique solution to 

(11.7) y = C (Dy) + C(-6y + Dy + f) 
we will have 

(11.8) y ± -► y. 

To see (11.8) note that ly i » 00 is bounded uniformly 
in i (from (11.6) and \i < 1) and if is defined as 

(11.9) g ± * sup ly-y I 

r >i 

Then {g^} form a nonincreasing sequence satisfying 

(11.10) g i+1 1 V g* + o(D 

1 x i-^co 

and (11.10) implies g i -*• 0 which, in turn, implies (11.8). 
The fact that y = y follows exactly as in Section 7 (see 
(7.27) ff) . Observe that since y is AP, By is also AP 
(from (11.1)). 

Now consider the operators defined by (8.6). The 
equation 

(11.11) Y h = c h (Dy h ) + C h (f) 
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can be solved exactly as (11.3) can. (For simplicity we 
suppose that B C^l! = II Cll so that we have the same restric- 
tion on £q-) 

Our first task is to show that y^ -*■ y. This will be a 
simple consequence of the strong convergence of to C. 

In fact if w^ = y - y^ we have 

(11.12) w h = [C-C h ](Dy) + C h (Dw h ) + [C-C h ] (f) 
whence 

(11.13) — ^Wj^oo + 

and this certainly implies II w. II > 0. 

n h^O 

Now of course this would not be useful unless we could 

prove 

(11.14) II w, II = OOi^) 

n 00 

This estimate can be shown provided we make certain 
assumptions on the interaction of D and the operator L 
defined in (8.11). We first observe that w^ satisfies 
the equation 

(11.15) w h = C h (Dw h ) + C h ([B-B h ]y) 

and using the analogue of (11.4) for the equation (11.11), 
together with the uniform boundedness of the operators , 
we see 

(11.16) = 0(» [ B - B h ]yS TO ) , 
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and it is sufficient to show 
< 11*17) I £B— B h ] y II = 0(h j ) * 

To obtain (11*17) we make an assumption on the operator 
L* We will call it Assumption I as we Will later have to make 
the same type of assumption for other operators. 

Assumption I: Suppose there exists a sequence of 

operators Lq = I, = i» such that 



spatial variable). 

Now Assumption I implies that if L^f is AP then L^y is AP , 
for & = 1 , . . . >n . First note that if z = C(f) then 


(11.19) = C(L £ f) . 

This is a consequence of (11. 18a, c) using the same argument 
as Was used for B (see (7.34) f f ) . Now if we define z £;i = L^y ± 
where the y^ are defined in (11.5); then it certainly follows 
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is AP for i and for £ « 0,..., 


from (11.18d) that z. . 

Ki / 1 

We will now show that the sequence {z„ . } converges 

J6 J 1 

To see this first set £ = 1. We then have 


n. 


as i-»- 


(11.20) 

Z, ... = C (Dz. . 
l;i+l l;i 

) + C(Ej(t) yi ) + C(f) 

and z. . 
J- / J- 

■+■ z 1 , where z 

is the solution to 

(11. 21) 

z 1 = CtDZj^) 

+ C(Ejy) + C(f) . 

In fact 

w l,i = Z 1 " z l,i 

satisfies the equation 

(11.21) 

w i,i.i = c(Dw 

l,i } + C(Ej(y-y.)) 

and 



(11.22) 

s: 

H- 

8 

11 

N 



follows from the same argument used in proving y is 
(see (11.9) ff ) . Finally by the closure of we must 
have = L^y. It is clear that a simple induction, using 
the same proof, will establish that 


(11.2 3) ^£^i — j_^oo > ' £ - l,...,n. 

In particular L (= L n )y -*- s ^ (11.17), hence (11.14) 

is valid. 

Finally we point out that if we assume a chain of 
operators, as in Assumption I, for the operator B, we will 
then have BDy AP. We will also have to make such an 
assumption for the operator L (see (10.16)), when applying 
the Lax-Wendroff scheme to (11.1), as will be discussed 
shortly . 



We will now apply to (11.1) the two methods of time 
discretization that have already been described. We are 
interested in the existence of AP sequence solutions, or 
equivalently exponential stability for the homogeneous 
difference equation, with k and h restricted according to 
relations (9.4) and (10.2). As one might expect this will 
be a simple consequence of Theorem 2, provided is 
sufficiently small. 

We consider first the linear multistep scheme as 

discussed in Section 9. Since the linear term C of the 

n 

difference scheme (see (2.11)) is a smooth function of its 

£+1 arguments, it follows that if k + k/h m is sufficiently 

small (i.e. UkB. fl is small) then the difference equation 
h 

for (11.1) can be written 

(H.24) » n+ l " C > n + k 0( Eo )w n + k f n 

where C* is the linear term corresponding to the unperturbed 
n 

version of (11.1), i.e. 
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of k and h (see (9.7) ff) it will follow from Theorem 2 

(see the remark following (9.16)) that in the same region 

the homogeneous version of (11.24) will be exponentially 

stable if is small enough. It then follows easily 

that the inhomogeneous equation (11.24) will have (L^,L^) 

~ n ~n 

(or (A ,A )) admissible, with an admissibility bound 
independent of k and h. 

We can thus obtain AP sequence solutions to (11.24) 
and, as shown in Section 9 (see (9.18) ff) , we will obtain 
convergence of 0(k p ) provided is bounded in h. 

However, as one can trivially see, if D, 6, and f are C p+ ^ 
then the same argument used to show y is (see (11.6) ff) 
can be repeated p+1 times and applied to the operators 
as well as C to show that y and y^ are C P+ ^ and in fact (11.7) 
can be generalized to (compare with (9.18)) 

(11.28) y 1 = C (Dy 1 ) + C < [ (D-<5 ) y ] 1 - (D - <S ) y 1 + f 1 ) , 

l ™ 1,..., pt 1 , 

with a similar formula for y^ and C^. This certainly shows 
that Uy^ll will be bounded in h for i = l,...,p+l and 
thus that we can obtain a complete extension of the theory 
of Section 9. We note in passing, that the same proof 
which led to (11.13) will show that y^ -+ y 1 . 

The extension of the Lax-Wendroff scheme will be just 
as simple. In fact (tinder the assumption that D, <5 and f 

3 

are C ) the difference equation for (11.1) will have the 
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same linear term as the equation for (11.25), (10.6) 
except for a perturbation 


(11.29) k [D + 


n 


I << B ch 


- <S. ) D 
n n 


D (B , 
n ch 


- S ) 
n 


+ D 


n 


+ °n)5 


as one can easily see by carrying out the expansion of 
Section 10 for the equation (11.1). It follows that if 
k/h" 1 - 0(1) (so that kllB^II is bounded). Theorem 2 is 
applicable, provided k and £q are sufficiently small. One 
can then obtain AP sequence solutions and carry out the 
convergence argument of Section 10 (see (10.17) ff) , 
with the only difficulty being the verification that 

(11.30) B (B 2 -B^)yll co = 0(h-^) 


and this will follow immediately if one assumes a chain as 
in Assumption I for the operator L introduced in (10.16 ) * 
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Appendix. Proof that Admissibility Implies 

an Exponential Dichotomy. 

Here we would like to give a simplified proof of the 
proposition that (A, A) admissibility for the inhomogeneous 
equation 

(A. 1) y = Ay + f 

implies an exponential dichotomy for the homogeneous equation 
(A. 2) y = Ay 

It is assumed of course that A(t) is AP although we will 
actually do the case. The proof is valid only if the 
underlying space E is finite dimensional and we can then 
regard A as a matrix and y and f as m vectors, where m is 
the dimension of E. The reader is referred to Section 5 
where an entirely similar proof for the difference equation 
case is given. 

The general proof is given in M & S (p. 344, Theorem 
10 3. A) , but this proof relies strongly on preceding material. 
A finite dimensional proof is given in Coppel (p. 134 ff) 
for an equation on [0,°°). His proof can be extended to an 
equation on the whole real axis, but the proof given here 
is simpler and more in keeping with the ideas of Massera 
and Schaffer. 

It is shown in M&S (Theorem 103. A) that (A, A) admissi- 
bility implies (L^,!^) admissibility and this will be our 
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starting point. We thus assume that for every f in L ot , (A.l) 
has a unique solution y in and the estimate 

(A. 3) II yB . < Kll fll 

will hold for a certain constant K. 

Wow let be the subspace of initial data which gives 
rise to solutions to (A. 2) that are bounded for t G [O, 00 )* 

Let S 2 be the analogous subspace giving rise to solutions 
that are bounded for t e (-“>,0]. We must have n {0} 

because any nontrivial/ bounded solution to (A. 2) would 
violate the uniqueness requirement of our definition of 
admissibility. Let Sy be any subspace complementary to 
so that 

(A. 4) E * S 1 © S 2 © S 3 

Let P^,P 2 ,P 3 be the associated projections, i.e. P^ is the 
projection onto S 1 along S 2 © S 3 , etc. We have 

(A. 5) I = P 1 + P 2 + *3 

P^P j —0, i = 1,2,3; j = 1,2,3? i ^ j. 

We will first show S 3 = {0}, i.e. P 3 = 0. This will 
follow from the variation of constants formula, just as in 
the difference equation case. If f has compact support and 
y(t) is the unique L m solution to (A.l) , then since y is a 
bounded solution to (A. 2) for large |t|, we must have 
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(A. 6) 


(a) 


(b) 


(c) 


(d) 


P 2 y 0 = 
P 3 y 0 = 
P l y 0 = 
P 3 y 0 = 


oo 





P 2 Y _1 (s) f { S ) 
P 3 Y _1 (s) f ( s ) 
P 1 Y~ 1 (s) f(s) 
P 3 y' 1 (s) f ( s ) 


ds 


ds 


ds 


ds 


Here Y(t) is the fundamental solution to (A. 2). 

Now for any vector z e E, if we set f = x^q jjY(t)z 
then (A.6b,d) yield immediately P 3 z = 0 and we can conclude 
P 3 = 0 , i.e. 

E = S ffi S, 

(A. 7) 1 

1 - P 1 + P 2 

Using (A.6a,c) we also see that for any f with compact 
support we have 


C X > 


(A. 8) 

y (t) 

= | G(t,s) f{s) ds 


where 






Y ( t ) P , Y _1 (s ) , 

t >_ S 

(A. 9) 

G(t,s) 

_ -L 

-Y(t)P 2 Y _1 (s) , 

S > t 


We will next use (A. 3) to show that for every fixed t/ 
G(t,s) is in (-oo f oo) as a function of s, and in fact 
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CO 


(A. 10) j ds ilG(t,s)H < K 

— OO 

where K is given in (A. 3). This is a simple consequence of 
using (A. 8) as an operator on L b (I), where I is any compact 
interval to conclude that 

(A. 11) j ds llGtt,s)« < K 

I 

and then letting I + (-<»,<») . 

We now will show that (A. 10) implies 

(A. 12) 0 G(t,s) 0 * 0(1) . 

To see this we simply observe that Y ^(s) satisfies the 
adjoint equation 

(A. 13) Z* = - Z A(s) 

where we have used to stand for d/ds. Hence fixing t 

and considering the region s £ t, we have 

(A. 14) G ( t , s ) = - G ( t / s ) A ( s ) 

Now since A is bounded it follows that G' (t,s) is in L^, as 
a function of s, on the interval (-°°,t3 with an norm 
which can be bounded uniformly in t. This of course means 
that G(t,s) must approach a limit as s -»• -«> and since G(t,s) 
is in (as a function of s) , this limit can only be zero. 
We thus have 

o 

(A- 15) G(t,s) = | dr G' (t,r) = 0(1) , s < t, 

— 03 
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and a similar argument establishes (A. 12) for s >_ t. 

We must now show that (A. 12) can be replaced by an 
exponentially decaying factor, and to do this we consider 
perturbations of (A. 1) and (A. 2); namely for small e > 0 
consider the equations 

(A. 16) y = [A-eI]y + f 

(A. 17) y = [A-eI]y 

If e is sufficiently small we will have (L ,L ) admissible 
for (A. 16). This is a simple consequence of the contracting 
mapping principle analogous to the argument of Section 4C (see 
(4.9) ff). The contracting mapping principle will also show that 
the admissibility bound for (A. 16) will be given by (see A. 3) 

(A. 18) K c - -■ K , < K , |ej < e n , 

e (1-eK) U 

if is sufficiently small. Thus the admissibility bound 

for (A. 16) can be taken independent of e, if e is small 
enough . 

Now the fundamental solution to (A. 18) is 
(A. 19) Y £ (t) = e~ £t Y ( t) 

and if P^, denote the projections replacing and # 
while G e (t,s) is the new Green's function, then the analysis 
which led to (A. 12), and which depended only on admissibility, 
is equally valid here and we can conclude 
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(A. 24) BP^ - P x « = 0<e) . 

To obtain (A. 22) we set H £ (t,s) = G e (t,s) — G(t»s). 
Now for any f with compact support we let y and y £ be the 
unique bounded solutions to (A.l) and (A. 1-6) respectively- 
Ifw = y - y we see that w is the unique bounded 
solution to 

(A. 25) w £ = Aw £ -e y £ , 
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and we can conclude from (A. 3) together with its analogue 
for (A. 16) that 

(A. 26) Hw £ II ^ « 0(e) II fl w 

Now w can also be written as 

^ 00 

(A. 27) w £ (t) * I ds H e (t,s) d (s) 

— 00 

and using the same argument that led to (A. 10) we see that 

oo 

(A. 28) | ds llH e (t,s)ll = 0(e) . 

— CO 

(Note that the bound in (A. 28) will be independent of t, 
although we will only use it for t = 0.) 

£ 

Now if s < t we see that as a function of s, H (t,s) 
satisfies the equation 

F ' £ £ 

(A. 29) H = - H A + e G (t,s) . 


This equation, together with (A. 20) and (A. 28), shows that 


(A. 30) | ds llH e (t, s) I = 0(e) 

— CO 

£ 

and since H (t,s) is in as a function of s we can write for 
s < t 

.31) H e (t,s) * 



and if we set t = 0, s = 0_ , (A. 31) yields (A. 24), which 
proves one part of the requirement of an exponential dichotomy 
and the other part follows on replacing -e by +e. 
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